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abstract  An  overview  of  the  parametrization  of  gravity  wave  drag  in  numerical  weather  prediction  and  climate 
simulation  models  is  presented.  The  focus  is  primarily  on  understanding  the  current  status  of  gravity  wave  drag 
parametrization  as  a  step  towards  the  new  parametrizations  that  will  be  needed  for  the  next  generation  of  atmos¬ 
pheric  models.  Both  the  early  history  and  latest  developments  in  the  field  are  discussed.  Parametrizations  devel¬ 
oped  specifically  for  orographic  and  convective  sources  of  gravity  waves  are  described  separately,  as  are  newer 
parametrizations  that  collectively  treat  a  spectrum  of  gravity  wave  motions.  The  differences  in  issues  in  and 
approaches  for  the  parametrization  of  the  lower  and  upper  atmospheres  are  highlighted.  Various  emerging  issues 
are  also  discussed,  such  as  explicitly  resolved  gravity  waves  and  gravity  wave  drag  in  models,  and  a  range  of 
unparametrized  gravity  wave  processes  that  may  need  attention  for  the  next  generation  of  gravity  wave  drag 
parametrizations  in  models. 

resume  [Traduit  par  la  redaction]  Le  present  article  donne  un  apergu  de  la  parametrisation  de  la  trainee  due 
aux  ondes  de  gravite  dans  les  modeles  numeriques  de  prevision  meteorologique  el  de  simulation  du  climat.  II  est 
axe  principalement  sur  la  comprehension  de  I’etat  actuel  de  la  parametrisation  de  la  trainee  due  aux  ondes  de 
gravite  en  vue  des  nouvelles  parametrisations  qui  seront  necessaires  pour  la  prochaine  generation  de  modeles  de 
l' atmosphere.  Les  debuts  ainsi  que  les  developpements  les  plus  recents  de  la  science  sont  examines.  Les 
parametrisations  developpees  expressement  pour  les  ondes  de  gravite  orographiques  et  les  ondes  de  gravite 
resultant  de  la  convection  sont  decrites,  de  meme  que  les  parametrisations  plus  recentes  visant  collectivement  une 
gamine  de  mouvements  des  ondes  de  gravite.  Les  differences  dans  les  preoccupations  et  les  methodes  propres  a 
la  parametrisation  de  T atmosphere  inferieure  et  de  T atmosphere  superieure  sont  soulignees.  Des  preoccupations 
nouvelles  sont  examinees,  notamment  les  ondes  de  gravite  et  la  trainee  due  aux  ondes  de  gravite  explicitement 
resolues  dans  les  modeles,  ainsi  qu’une  gamine  de  processus  associes  aux  ondes  de  gravite  qui  ne  sont  pas 
parametrises  mais  dont  ilfaudra  eventuellement  tenir  compte  dans  la  prochaine  generation  de  parametrisations 
de  la  trainee  due  aux  ondes  de  gravite  dans  les  modeles. 


1  Introduction 

There  have  been  considerable  advances  in  the  numerical  pre¬ 
diction  of  weather  ever  since  Lewis  Fry  Richardson  attempt¬ 
ed  the  first  weather  prediction  during  World  War  I.  With  the 
advance  of  human  knowledge  and  the  advent  of  state-of-the- 
art  computing  technologies,  it  is  becoming  more  feasible  to 
tackle  the  formidable  problem  of  accurate  numerical  weather 
and  climate  prediction. 

Although  the  efforts  to  solve  the  problem  of  numerical 
weather  prediction  (NWP)  are  severely  hindered  by  limita¬ 


tions  in  physical  theories,  mathematical  techniques  and  com¬ 
putational  expense,  significant  developments  have  been 
achieved  in  all  the  areas  of  the  atmospheric  sciences.  Many 
concurrent  efforts  are  being  made  to  represent  the  physical 
processes  of  the  atmosphere  more  accurately  and  more  effi¬ 
ciently,  along  with  parallel  efforts  to  solve  the  system  of 
equations  governing  the  atmosphere  more  accurately  and 
more  efficiently  in  NWP  and  general  circulation  models 
(GCMs).  Contemporary  NWP  and  climate  models  are  now 
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successful  in  simulating  major  observed  mean  features  of  the 
atmosphere.  However,  even  the  most  up-to-date  models  often 
suffer  from  some  unrealistic  aspects  of  their  simulated  atmos¬ 
pheres.  Notable  examples  are  the  “cold-pole”  problems  asso¬ 
ciated  with  an  unrealistically  strong  polar  night  jet  in  the 
stratosphere  (e.g..  Shepherd,  2000)  which  is  closely  linked  to 
excessively  zonal  and  strong  surface  westerlies  (also  called 
the  “westerly  bias”),  as  first  noted  in  the  northern  hemisphere 
winter  (e.g..  Palmer  et  al.,  1986). 

Aside  from  the  theoretical  deficiencies,  even  the  most  pow¬ 
erful  available  computing  architectures  still  cannot  run  typi¬ 
cal  NWP  or  climate  models  fast  enough  to  resolve  all  relevant 
scales  of  atmospheric  motion.  At  present,  global  models 
must,  in  practice,  be  run  with  horizontal  resolutions  that 
cannot  typically  resolve  atmospheric  phenomena  shorter 
than  -10-100  km  or  greater  for  weather  prediction  and 
-100-1000  km  or  greater  for  climate  prediction.  Many  atmos¬ 
pheric  processes  have  shorter  horizontal  scales  than  these, 
and  some  of  these  “subgrid-scale”  processes  interact  with  and 
affect  the  larger-scale  atmosphere  in  important  ways.  Since 
they  cannot  be  resolved,  large-scale  models  must  resort  to 
“parametrizations”  that  capture  the  salient  effects  on  the 
resolved  atmosphere. 

Atmospheric  gravity  waves  are  one  such  unresolved 
process.  These  waves  are  generated  by  lower  atmospheric 
sources,  e.g.,  flow  over  irregularities  at  the  Earth’s  surface 
such  as  mountains  and  valleys,  uneven  distribution  of  diabat- 
ic  heat  sources  associated  with  convective  systems,  and  high¬ 
ly  dynamic  atmospheric  processes  such  as  jet  streams  and 
fronts.  The  dissipation  of  these  waves  produces  synoptic- 
scale  body  forces  on  the  atmospheric  flow,  known  as  “gravi¬ 
ty-wave  drag”  (GWD),  which  affects  both  short-term 
evolution  of  weather  systems  and  long-term  climate. 
However,  the  spatial  scales  of  these  waves  (in  the  range  of 
-5-500  km  horizontally)  are  too  short  to  be  fully  captured  in 
models,  and  so  GWD  must  be  parametrized.  It  is  now  gener¬ 
ally  agreed  that  the  westerly  bias  in  models  can  be  alleviated, 
to  some  extent,  with  a  suitable  parametrization  of  GWD  (e.g., 
Boville,  1995).  While  the  adequacy  of  some  aspects  of  the 
formulations  (and  even  the  necessity  of  such  a  parametriza¬ 
tion  for  westerly  bias)  is  still  in  debate  (Pawson  et  al.,  1998; 
Newman  et  al.,  2001),  the  role  of  GWD  in  driving  the  global 
middle  atmosphere  circulation  and  thus  global  mean 
wind/temperature  structure  is  well  established.  Thus,  GWD 
parametrizations  are  now  critical  components  of  virtually  all 
large-scale  atmospheric  models. 

The  parametrization  of  the  effects  of  atmospheric  gravity 
waves  has  evolved  into  branches  with  their  own  (interrelated) 
backgrounds,  histories  and  communities.  This  paper  attempts 
to  provide  a  global  overview  of  the  past,  present  and  foresee¬ 
able  future  in  the  parametrization  of  gravity  wave  effects  in 
NWP  and  climate  models,  with  the  emphasis  on  GWD. 
Throughout  this  paper  we  will  continue  to  use  the  common 
term  “gravity  wave  drag”  to  describe  the  forces  exerted  on  the 
atmosphere  by  gravity  wave  dissipation,  despite  its  shortcom¬ 
ings:  “drag”  implies  a  deceleration  of  the  flow,  whereas 


the  body  forces  on  the  flow  produced  by  gravity  wave  dissi¬ 
pation  can  either  accelerate  or  decelerate  (retard  or  reverse) 
atmospheric  winds.  This  broader  definition  of  GWD  should 
be  kept  in  mind  throughout  this  article  (Appendix  Ad3  has 
further  details). 

Section  2  discusses  the  early  history  of  and  motivations 
behind  GWD  parametrization.  Section  3  reviews  parame¬ 
trizations  of  gravity  waves  from  two  important  sources: 
mountains  and  convection.  Section  4  discusses  newer  param¬ 
etrizations  originally  designed  for  middle  atmospheric  models 
that  are  now  being  applied  more  widely.  Methods  of  evaluat¬ 
ing  these  parametrization  schemes  are  discussed  in  Section  5. 
Since  some  gravity  waves  are  being  resolved  in  new  high-res- 
olution  models,  we  review  some  of  these  findings  in 
Section  6.  Section  7  discusses  a  number  of  specific  issues  that 
may  need  to  be  considered  in  future  development  work.  The 
overall  situation  is  briefly  summarized  in  Section  8. 

Given  that  this  is  an  overview  paper,  many  details  will  nec¬ 
essarily  be  omitted.  Appendix  A  provides  a  brief  tour  through 
common  nomenclature,  concepts  and  dynamical  theories  of 
gravity  waves  relevant  to  GWD  parametrization.  Those  read¬ 
ers  who  are  new  to  this  topic  may  wish  to  read  this  Appendix 
before  proceeding  to  the  next  section.  In  addition,  references 
are  provided  to  the  relevant  parts  of  Appendix  A  in  the  fol¬ 
lowing  sections,  to  aid  readers  who  may  need  theoretical 
background  on  certain  topics.  Readers  needing  more  detailed 
background  on  atmospheric  gravity  waves  can  find  it  in  vari¬ 
ous  texts  (e.g.,  Gossard  and  Hooke,  1975;  Gill,  1982; 
Andrews  et  al.,  1987;  Baines,  1995)  and  review  articles 
(Smith,  1979,  2001;  Fritts,  1984;  Hamilton,  1998;  Fritts  and 
Alexander,  2003).  Additionally,  the  basic  ways  in  which 
GWD  drives  global  atmospheric  circulation  features  have 
been  described  quite  well  in  a  number  of  review  articles  (e.g., 
Andrews,  1987;  Holton  and  Alexander,  2000;  McIntyre, 
2001).  Subsequent  sections  and  Appendix  A  also  provide 
citations  to  key  review  papers  on  specific  subject  matter  in 
greater  detail.  Finally,  Appendix  B  lists  the  acronyms  and 
symbols  used  in  this  paper. 

2  Birth  of  gravity  wave  drag  parametrization 

a  Misrepresented  Damping  Mechanisms  in  Models 
The  earliest  large-scale  models  with  relatively  low  resolutions 
and  low  model  tops  seemed  to  simulate  observed  features  of 
the  subtropical  jet  quite  well,  at  least  for  the  northern  hemi¬ 
sphere  winter  (e.g.,  Kasahara  et  al.,  1973).  With  increased 
spatial  resolution,  however,  the  simulated  jet  became  exces¬ 
sively  strong.  Studies  initially  investigated  whether  underes¬ 
timated  surface  friction  drag  contributed  mostly  to  this  wind 
bias  (Swinbank,  1985).  If  friction  drag  is,  say,  underestimat¬ 
ed,  surface  drag  due  to  grid-scale  mountains  can  take  over 
what  should  have  been  represented  by  friction  drag  and  in 
turn  be  overestimated  in  a  model  in  order  to  maintain  the 
global  angular  momentum  balance.  Studies  soon  concluded 
that  underestimated  friction  drag  was  not  the  main  cause  of 
the  bias,  mainly  because  friction  drag  is  not  easily  enhanced 
under  the  stable  atmospheric  conditions  of  winter  for  which 
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Mid-latitude  December  (Northern  Winter)  Mid-latitude  July  (Northern  Summer) 


zonal  winds  U(z)  zonal  winds  U(z) 


gravity  wave  phase  velocity  c  gravity  wave  phase  velocity  c 

(m  s_1  eastward)  (m  s_l  eastward) 

Fig.  1  Typical  mid-latitude  zonal  winds  U(z )  during  northern  (a)  winter  and  (b)  summer.  Black  curve  shows  observed  winds,  grey  curve  shows  model  “radia¬ 
tive”  winds  that  result  without  a  wave  drag  parametrization.  Sources  of  gravity  waves  with  various  phase  speeds  c  are  also  depicted,  with  the  source 
and  wave  breaking  symbols  similar  to  those  defined  in  Fig.  10.  On  these  plots,  waves  ascend  vertically  upwards  since  c  remains  constant,  until  they 
break  or  reach  a  critical  level  c  =  U(zc )•  (Based  on  a  presentation  first  used  by  Lindzen,  1981) 


the  bias  was  most  dominant.  Moreover,  global  increases  in 
surface  friction  drag  can  have  deleterious  effects  elsewhere: 
e.g.,  excessively  weak  southern  hemispheric  circumpolar 
flow  (Boer  and  Lazare,  1988).  (For  a  discussion  on  various 
drag  processes  that  arise  at  the  surface,  see  Appendix  Aa.) 

If  the  atmosphere  were  in  radiative  equilibrium,  there 
would  exist  a  polar  night  jet  in  the  middle  atmosphere  that  is 
much  stronger  and  extends  higher  than  is  observed  (see 
Fig.  1;  see  also  Fig.  6  of  Hamilton,  1996).  On  the  other  hand, 
“Rayleigh  friction”  (e.g.,  Boville,  1986)  and  “Newtonian 
cooling”  (Dickinson,  1973),  which  exponentially  damp  winds 
and  temperatures,  respectively,  to  some  reference  state  as  a 
function  of  time  and  are  often  used  to  suppress  spurious 
reflection  of  waves  from  model  tops,  were  found  to  be  quite 
effective  in  reducing  the  magnitude  of  the  polar  night  jet  in 
models  (Leovy,  1964;  Schoeberl  and  Strobel,  1978;  Holton 
and  Wehrbein,  1980).  The  success  here  is  qualified:  an  impor¬ 
tant  drawback  of  Rayleigh  damping  is  that  it  always  drags  the 
winds  back  toward  zero  (or  climatology),  whereas  observed 
winds  go  to  zero  then  generally  change  direction  (see  Fig.  1). 
Moreover,  this  kind  of  ad-hoc  damping  usually  requires  a 
thick  vertical  model  domain  and  may  produce  excessive 
damping  accompanied  by  unrealistic  changes  and  reductions 
in  atmospheric  variability  (e.g.,  Shepherd  et  ah,  1996; 
Lawrence,  1997a;  Kim  et  ah,  1998).  Nonetheless,  the  ability 
of  such  ad-hoc  damping  over  a  thick  stratospheric  domain  to 
alleviate  the  problem  of  excessively  strong  stratospheric  jets 
suggested  that  some  kind  of  drag  mechanism  was  missing 
from  the  models. 

Later  studies  eventually  concluded  that  the  wind  bias  was 
mainly  due  to  the  lack  of  an  explicit  simulation  of  “drag”  gen¬ 
erated  by  breaking  subgrid-scale  gravity  waves  (Houghton, 
1978;  Lindzen,  1981;  Matsuno,  1982;  Holton,  1982,  1983; 


details  are  discussed  in  Section  3).  Earlier  lower  top  models 
enjoyed  apparently  successful  simulation  of  the  jet  because 
the  resolution  was  coarse  enough  for  the  underestimated 
meridional  eddy  momentum  transport  by  planetary  waves  to 
balance  approximately  the  underestimated  vertical  momen¬ 
tum  transport,  as  illustrated  in  Fig.  2.  When  the  horizontal  res¬ 
olution  of  these  models  increased,  the  meridional  momentum 
transport  was  more  accurately  resolved  and  thus  increased, 
and  it  could  no  longer  be  balanced  by  the  still  underestimated 
vertical  momentum  transport.  Without  a  parametrization  of 
subgrid-scale  GWD,  the  mid-latitude  westerly  (easterly) 
winds  in  winter  (summer)  became  excessively  strong  due  to 
unresolved  momentum  transport  to  the  ground  that  should 
balance  the  resolved  meridional  momentum  transport  (see 
Section  7  of  Palmer  et  al.,  1986  and  Section  5  of  Boer  and 
Lazare,  1988). 

b  Enhanced  Orography  with  Effects  Resembling  Gravity 
Wave  Drag 

Although  there  were  several  pioneering  theoretical  and  ana¬ 
lytical  studies  on  the  treatment  and  effects  of  GWD  on  the 
large-scale  background  flow  (e.g..  Sawyer,  1959;  Lindzen, 
1981;  Holton,  1982),  the  first  practical  attention  to  this  sub¬ 
ject  came  with  the  introduction  of  “envelope  orography” 
(e.g.,  Wallace  et  al.,  1983).  Envelope  orography  is  a  type  of 
grid-scale  orography  in  which  mountain  heights  are  elevated 
proportionally  to  the  standard  deviation  of  the  subgrid-scale 
orographic  elevations  within  each  grid  box.  Another  similar 
type  of  orography  is  “silhouette  orography”  in  which  model 
orography  is  given  by  the  average  of  the  mountain  peaks 
within  each  grid  box  (Mesinger  and  Janjic,  1986):  Mesinger 
and  Collins  (1986)  and  Lott  and  Miller  (1997b)  review  vari¬ 
ous  model  representations  of  grid-scale  orography.  By  using 
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Fig.  2  Simplified  schematic  of  the  angular  momentum  transport  between  the  atmosphere  and  the  Earth’s  surface  in  atmospheric  models  (northern  hemisphere). 
The  angular  momentum  is  transported  vertically  from  westerlies  to  the  ground  over  a  mountainous  area  in  middle  latitudes  and  from  the  ground  to  east¬ 
erlies  in  low  latitudes.  To  balance  the  loss  of  momentum  in  the  westerlies  region  and  the  gain  in  the  easterlies  region,  meridional  momentum  transport 
arises  to  maintain  mass  continuity. 


these  elevated  forms  of  orography  in  large-scale  models,  it 
was  shown  that  the  biases  in  the  zonal  mean  wind  and  tem¬ 
perature  were  reduced  through  enhanced  “mountain  drag” 
(see  Appendix  Aa),  i.e.,  enhanced  generation  of  planetary 
wave  activity  and  associated  enhanced  eddy  flux  divergences 
that  impact  mean  fields  (e.g.,  Wallace  et  al.,  1983;  Palmer  and 
Mansfield,  1986;  Tibaldi,  1986;  Iwasaki  and  Sumi,  1986). 

The  improvements  gained  by  enhancing  the  orography 
came  via  greater  generation  of  Rossby  waves,  and  not  from 
parametrized  effects  of  subgrid-scale  gravity  waves  (even 
though  information  on  the  subgrid-scale  orography  is  used  to 
create  the  grid-scale  orography).  While  this  is  of  course  true, 
enhanced  orography  can  be  viewed  in  a  different  context  in 
terms  of  the  spectrum  of  near-surface  gravity  waves  associat¬ 
ed  with  subgrid-scale  orography  (Kim,  1996).  Enhancing 
orography  is,  in  a  sense,  equivalent  to  considering  stagnant 
flow  formed  over  mountainous  terrain,  which  is  relatively 
deep/shallow  over  valleys/ridges  (Fig.  3).  This  stagnant  flow 
effectively  acts  as  a  barrier  or  “envelope”  to  the  grid-scale 
flow.  The  streamlines  of  the  terrain  and  the  elevated  surfaces 
due  to  stagnant  flow  resemble  those  of  “external  gravity 
waves”  (see  Appendix  Acl;  Fig.  11)  that  are  evanescent,  i.e., 
decay  in  the  vertical. 

Envelope  orography  has  been  investigated  in  relation  to  its 
impact  on  synoptic-scale  blocking  (e.g.,  Mullen,  1994).  Some 
recent  models  use  explicit  parametrizations  of  lift  forces  due 
to  blocking  near  the  surface  complementing  the  work  by 
GWD  (Fott,  1999).  Caution  is,  however,  needed  to  distin¬ 
guish  between  the  drag-decreasing  effect  of  upstream  flow¬ 
blocking  and  the  drag-enhancing  effect  of  downstream 


low-level  wave -breaking  (Kim  and  Arakawa,  1995;  Kim  et  al. 
1998;  see  also  Fig.  4).  For  example.  Smith  et  al.  (2002) 
recently  showed  how  a  low-level  stagnant  layer  reduced  the 
upward  generation  of  (internal)  mountain  waves  (see  Section 
7e2  for  further  discussion). 

One  of  the  drawbacks  of  enhanced  orography  is  that  these 
stagnant  flow  effects  are  parametrized  in  an  invariant  way, 
whereas  in  reality  they  are  time  (i.e.,  flow)  dependent  (Fott 
and  Miller,  1997b).  Further,  envelope  orography  may  inter¬ 
fere  with  data  assimilation  by  necessitating  rejection  of  low- 
level  atmospheric  data,  and  may  also  generate  excessive 
precipitation  in  models  (Fott  and  Miller,  1997a).  There  has 
been  some  debate  on  the  relationship  between  and  the  mech¬ 
anisms  behind  enhanced  orography  and  GWD  (e.g..  Section  1 
of  Kim,  1996).  With  successive  increases  in  horizontal  reso¬ 
lutions,  however,  newer  models  have  tended  to  use  more  real¬ 
istic  “mean”  grid-scale  orography  instead  of  enhanced 
orography.  It  seems  that  the  earlier  beneficial  impact  of 
enhanced  orography  was  partially  due  to  inadequate  horizon¬ 
tal  resolution. 

3  Evolution  of  gravity  wave  drag  parametrization 

Studies  subsequent  to  those  discussed  in  the  previous  section 
soon  evolved  into  explicit  consideration  of  the  breaking  of 
vertically  propagating  “internal”  gravity  waves  (for  an  intro¬ 
duction  to  wave  breaking  and  GWD,  see  Appendix  Ad).  First 
parametrizations  were  for  gravity  waves  generated  by  flow 
over  subgrid-scale  orography  (also  called  mountain  waves), 
followed  by  parametrizations  of  horizontally  propagating,  but 
vertically  “trapped"  lee  waves  (see  Appendix  Abl  and  Ac  for 
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Fig.  3  Schematic  figure  showing  the  enhancement  of  large-scale  orography  due  to  the  formation  of  stagnant  flow  over  subgrid-scale  orography.  Large-scale 
flow  encounters  the  effective  height  of  the  orography  instead  of  the  average  height  since  it  does  not  flow  into  stagnant  regions  at  the  base  of  valley 
regions.  Envelope  orography  is  an  example  of  such  enhanced  orography.  The  streamline  patterns  resemble  those  of  external  gravity  waves  (i.e.,  that  are 
vertically  evanescent).  (Taken  from  lecture  notes  of  Prof.  Akio  Arakawa.) 
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Grid  Interval  of  Large-Scale  Model  Momentum  Flux 

Fig.  4  (Left)  Schematic  figure  showing  low-level  wave-breaking  and  stagnation  points  generated  in  a  model  by  resonant  amplification  of  mountain  waves 
between  the  surface  and  a  wave-induced  critical  layer  near  the  peak  of  the  mountains  downstream.  Contours  and  shades  represent  isentropes  and 
decreased  wind  regions  respectively.  Also  shown  are  a  region  of  upper-level  wave  breaking  due  to  the  breakdown  of  vertically  propagating  mountain 
waves,  and  a  region  of  flow-blocking  in  the  upstream.  (Right)  A  typical  profile  of  horizontally  averaged  vertical  momentum  flux  corresponding  to  the 
upstream  and  downstream  regions  of  the  schematic  figure  on  the  left.  (Schematically  drawn  based  on  Kim  and  Arakawa,  1995) 
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background  on  these  topics).  With  the  tops  of  the  models  now 
extending  upward,  a  new  need  has  emerged  to  parametrize 
drag  at  upper  levels  due  to  dissipation  of  “non- stationary” 
gravity  waves  (see  Appendices  Ab2  and  Ac);  issues  that  are 
discussed  specifically  in  Section  4.  In  tropical  regions,  gravi¬ 
ty  waves  and  tropical  planetary  waves  (e.g.,  Kelvin  waves) 
generated  by  deep  convection  are  considered  the  major  source 
of  non- stationary  wave  drag  that  drives  atmospheric  variabil¬ 
ity  in  the  stratosphere  and  above.  Since  orography  and  con¬ 
vection  are  currently  the  two  most  parametrized  gravity  wave 
sources,  we  discuss  the  evolution  of  parametrized  GWD  from 
both  of  these  sources  in  the  remainder  of  this  section. 

a  Orographic  Gravity  Wave  Drag  Parametrization 
In  the  presence  of  favourable  background  winds  and  static 
stability,  gravity  waves  are  generated  by  flow  over  mountain¬ 
ous  terrain.  The  literature  on  “mountain  waves”  is  extensive 
and  comprehensive  reviews  of  the  relevant  theory  are  avail¬ 
able  (e.g..  Smith,  1979,  2001;  Baines,  1995;  Wurtele  et  al., 
1996).  Modelling  of  mountain  waves  goes  back  many  years 
and  mathematical  theories  have  evolved  that  clarify  many 
aspects  of  wave  generation  and  evolution.  This  work  greatly 
aided  rapid  development  of  parametrization  schemes. 

The  first  generation  GWD  parametrization  schemes  for 
large-scale  models  with  relatively  low  model  tops  were  devel¬ 
oped  for  mountain  waves  (e.g.,  Boer  et  al.,  1984;  Palmer 
et  al.,  1986;  McFarlane,  1987).  The  first  formulations  were 
single-wave  parametrizations  based  on  two-dimensional  (2D) 
linear  single-wave  nonrotating  stationary  hydrostatic  gravity- 
wave  theory,  utilizing  (except  for  Boer  et  al.  1984)  the  “satu¬ 
ration  hypothesis”  (Lindzen,  1981):  see  Appendix  Ad  and 
Section  4b  for  further  details.  One  of  the  main  tasks  of  these 
schemes  was  to  separate  the  stratospheric  polar  night  jet  from 
the  tropospheric  subtropical  jet  by  reducing  the  overall  mag¬ 
nitude  of  the  jets  and  creating  stronger  easterly  wind  shear  in 
the  upper  troposphere.  A  typical  impact  of  these  schemes  is 
very  large  mid-latitude  lower  stratospheric  drag  that  has  a 
direct  impact  on  the  stratospheric  jet  and  an  indirect  impact 
on  the  surface  westerlies  through  the  secondary  circulation 
induced  by  the  stratospheric  drag  (Fig.  5).  This  indirect 
impact  was  effective  in  warming/cooling  the  polar/tropical 
stratosphere  (i.e.,  reducing  cold-pole  problems)  and  decreas¬ 
ing  the  surface  westerlies  associated  with  increased  zonal 
asymmetry  in  surface  pressure  fields  (i.e.,  reducing  the  west¬ 
erly  bias). 

At  about  the  same  time,  within  the  mesoscale  modelling 
community  there  were  active  discussions  on  the  mechanisms 
responsible  for  the  so-called  “severe  downslope  windstorms” 
found  downstream  of  major  mountains:  e.g.,  Boulder  wind¬ 
storms  downstream  of  the  Rockies  (Klemp  and  Lilly,  1978; 
Peltier  and  Clark,  1979;  Durran  and  Klemp,  1987).  An  illus¬ 
tration  of  flow  over  a  relatively  high,  isolated  mountain  based 
on  a  2D  mesoscale  model  simulation  of  this  situation  is 
shown  in  Fig.  4.  “Wave  resonance”  occurs  through  formation 
of  a  self-induced  mountain  wave  critical  level  associated  with 
mountain  wave  breaking  (see  Appendices  Ac3  and  Ad4)  that 


vertically  traps  and  then  amplifies  wave  energy.  Although 
there  has  been  some  debate  about  underlying  physical  mech¬ 
anisms  (based  on  linear  hydrostatic,  non-linear  non-hydrosta- 
tic  or  internal  hydraulic  theories),  it  is  generally  agreed  that 
the  GWD  associated  with  downslope  windstorms  can  be  very 
large  (e.g.,  Peltier  and  Clark,  1979;  Bacmeister  and 
Pierrehumbert,  1988).  While  these  events  are  intermittent,  the 
estimated  drag  each  event  produces  is  often  larger  than  that 
produced  by  stratospheric  mountain  wave  breaking.  Similar 
resonance  breaking  and  drag  may  also  be  important  in  the 
boundary  layer  (Nappo  and  Chimonas,  1992).  Thus  this 
process  is  important  for  the  large-scale  atmospheric  momen¬ 
tum  budget. 

Subsequent  development  of  orographic  GWD  parametriza¬ 
tion  was  motivated  by  findings  from  these  mountain  wave 
simulations  and  observations,  as  well  as  the  need,  in  models, 
for  increased  low-level  drag  (see  Section  7d).  Some  schemes 
began  to  represent  these  effects  through  a  simple  lower  tro¬ 
pospheric  enhancement  in  GWD  (e.g.,  Pierrehumbert,  1986). 
The  effects  of  “linearly  trapped"  non-hydrostatic  lee  waves 
downstream  of  mountains  were  parametrized  in  an  ad-hoc 
manner  by  Iwasaki  et  al.  (1989),  resulting  in  improved  fore¬ 
casts.  The  effects  of  “non-linearly  trapped"  lee  waves  in  the 
lower  troposphere  due  to  low-level  wave-breaking  were  sys¬ 
tematically  parametrized  to  generate  selective  enhancements 
(only  in  downstream  regions  with  strong  non-linearity)  of  low- 
level  drag  (Kim  and  Arakawa,  1995)  and  implemented  as  a 
parametrization  in  both  a  climate  simulation  model  (Kim, 
1996)  and  a  weather  forecast  model  (Alpert  et  al.,  1996),  with 
improvements  in  each  case.  This  selective  enhancement  dis¬ 
tinguished  between  the  low-level  wave -breaking  region  down¬ 
stream  and  the  flow-blocking  region  upstream  that  are 
associated  with  stronger  and  weaker  vertical  divergence  of 
horizontal  momentum  flux,  respectively  (Fig.  4).  The  latest 
mountain  wave  parametrization  schemes  have  devoted  much 
effort  to  developing  and  improving  these  and  other  lower-level 
drag  and  orographic  specifications  (e.g.,  Kim  and  Arakawa, 
1995;  Lott  and  Miller,  1997a;  Gregory  et  al.,  1998;  Scinocca 
and  McFarlane,  2000),  which  led  to  improvements  in  the  glob¬ 
al  models’  overall  forecast  skill  and  simulation  of  mean  sea 
level  pressures  (e.g.,  Milton  and  Wilson,  1996;  Kim,  1996). 

Various  versions  of  first-generation  orographic  GWD 
parametrization  schemes  are  now  routinely  implemented  in 
large-scale  models  for  both  climate  simulations  and  weather 
forecasts  (e.g.,  Boer  et  al.,  1984;  Palmer  et  al.,  1986; 
McFarlane,  1987;  McFarlane  et  al.,  1987;  Helfand  et  al.,  1987; 
Rind  et  al.,  1988;  Miller  et  al.,  1989;  Iwasaki  et  al.,  1989; 
Broccoli  and  Manabe,  1992;  Kim,  1996;  Milton  and  Wilson, 
1996;  Chun  et  al.,  1996;  Zhou  et  al.,  1996;  Lott  and  Miller, 
1997a;  Kiehl  et  al.,  1998;  Gregory  et  al.,  1998;  Lott,  1999; 
Scinocca  and  McFarlane,  2000).  Intercomparisons  of  various 
parametrization  schemes  in  models  are  available  in  the  litera¬ 
ture  (e.g.,  Boer  et  al.,  1992;  Gates,  1992;  Kim  and  Arakawa, 
1995;  Pawson  et  al.,  2000).  Orographic  GWD  parametrization 
is  in  a  state  of  continual  development.  A  range  of  issues  for  its 
future  development  is  set  out  in  Section  7. 
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Fig.  5  An  impact  of  a  stratospheric  momentum  sink  (with  its  maximum  at  45°N  and  limits  in  higher  and  lower  latitudes)  on  the  zonally  averaged  large-scale 
circulation  for  the  northern  winter.  The  stratospheric  drag  in  the  mid-latitude  is  deflected  to  the  North  Pole  due  to  Coriolis  force,  down  to  the  tropos¬ 
phere,  to  the  middle  latitudes,  and  back  to  the  stratosphere  in  the  low  latitudes,  while  inducing  adiabatic  warming/cooling  in  the  high/low  latitudes  and 
deceleration  of  the  winds  throughout  the  atmosphere  in  middle  latitudes.  (Schematically  constructed  based  on  Palmer  et  al.,  1986) 


b  Convective  Gravity  Wave  Drag  Parametrization 
While  more  attention  has  been  paid  to  the  parametrization  of 
orographic  gravity  waves  to  alleviate  the  wind  and  tempera¬ 
ture  biases  in  the  mid-latitude  northern  winter  hemisphere, 
significant  interest  has  recently  emerged  in  the  parametriza¬ 
tion  of  gravity  waves  generated  by  convective  systems.  The 
importance  of  convectively-generated  tropical  waves  in  dri¬ 
ving  the  equatorial  stratospheric  semi-annual  oscillation 
(SAO)  and  quasi-biennial  oscillation  (QBO)  has  been  appre¬ 
ciated  for  many  years,  but  the  relative  importance  of  the 
planetary-scale,  intermediate-scale  and  small-scale  (subgrid- 
scale)  waves  to  the  equatorial  momentum  budget  has  been 
debated  for  many  years  (e.g.,  Lindzen  and  Holton,  1968; 
Holton  and  Lindzen,  1972;  Hitchman  and  Leovy,  1988; 
Boville  and  Randel,  1992;  Bergman  and  Salby,  1994; 


Dunkerton,  1997;  Garcia,  2000).  Global  models  can  be  used 
to  test  this  to  some  extent,  since  planetary-scale  waves  can  be 
explicitly  simulated  while  GWD,  due  to  smaller-scale  waves, 
can  be  separately  parametrized. 

Simulations  of  the  QBO  and/or  SAO  in  three-dimensional 
(3D)  global  models  have  been  attempted  in  various  ways, 
using,  for  example;  a  mechanistic  model  that  either  specifies 
large  arbitrary  wave  amplitudes  (e.g.,  Takahashi  and  Boville, 
1992)  or  parametrizes  GWD  (e.g.,  Lawrence,  2001);  a  sim¬ 
plified  high-resolution  GCM  (Horinouchi  and  Yoden,  1998); 
a  full  GCM  with  either  a  GWD  parametrization  (e.g.,  Jackson 
and  Gray,  1994;  Scaife  et  al.,  2000),  increased  vertical  reso¬ 
lution  (Hamilton  et  al.,  1999,  2001),  or  an  additional  stratos¬ 
pheric  dry  convective  adjustment  process  (Takahashi,  1996). 
This  recent  work  with  full  GCMs  suggests  that  higher  fre- 
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Fig.  6  Gravity  waves  generated  by  penetrative  convection  in  the  presence  of  background  wind.  (Based  on  Hooke,  1986;  by  permission  of  the  American 
Meteorological  Society).  Similar  wave  features  are  seen  in  the  latest  numerical  model  simulations  (e.g.,  Holton  and  Alexander,  1999). 


quency  tropical  gravity  waves  generated  by  convection  may 
need  to  be  represented  better  in  large-scale  models  to  improve 
simulations  of  these  features  (see  Section  6a).  Thus,  parame- 
trizations  of  convectively-generated  GWD  can  be  considered 
a  means  of  alleviating  these  deficiencies,  allowing  better  sim¬ 
ulations  of  important  tropical  circulation  features  such  as  the 
QBO  in  the  lower  stratosphere  and  the  SAO  in  the  upper 
stratosphere  and  mesosphere. 

Compared  with  orographic  gravity  waves,  it  has  proven 
more  difficult  to  model  the  way  in  which  gravity  waves  are 
generated  by  various  convective  sources;  the  simplest  situa¬ 
tion  is  depicted  in  Fig.  6.  There  are  currently  three  proposed 
generation  mechanisms  in  the  literature.  The  first  considers 
convective  clouds  in  terms  of  their  “diabatic”  or  “thermal” 
forcing.  In  this  situation,  thermal  forcing  terms  within  the 
convective  clouds  radiate  gravity  waves  into  the  stably  strati¬ 
fied  atmosphere  above  the  clouds  (Salby  and  Garcia,  1987; 
Bergman  and  Salby,  1994;  Garcia,  2000).  It  has  been  argued 
that  the  vertical  wavelengths  of  convectively  generated  grav¬ 
ity  waves  are  approximately  twice  the  heating  depth  and  thus 
intrinsic  phase  speeds  of  these  waves  are  determined  mainly 
by  the  vertical  depth  of  the  heating  region  (Bergman  and 
Salby,  1994;  Alexander  et  al.,  1995;  Garcia,  2000).  In  con¬ 
trast,  Holton  et  al.  (2002)  showed  that  the  vertical  wavelength 
in  their  study  depends  on  a  non-dimensional  parameter  relat¬ 
ed  to  the  frequency,  horizontal  and  vertical  scales  of  the  forc¬ 
ing.  Detailed  mesoscale  model  simulations  by  Song  et  al. 


(unpublished  manuscript,  2003)  suggest  a  dependence  on  the 
basic-state  wind  relative  to  the  rearward  propagating  convec¬ 
tive  cells.  A  series  of  analytical  studies  of  thermally-induced 
internal  gravity  waves  related  to  the  thermal  forcing  mecha¬ 
nism  are  available  in  the  literature  (e.g..  Smith  and  Lin,  1982; 
Lin  and  Smith,  1986;  Lin,  1987;  Bretherton,  1988;  Chun, 
1995;  Chun  and  Baik,  1998;  Baik  et  al.,  1999a,b). 

The  second  mechanism  assumes  that  convective  heating  in 
clouds  produces  “obstacles”  to  the  flow  in  the  form  of  verti¬ 
cal  bulges  of  isentropic  layers  above  the  convective  region. 
Gravity  waves  are  generated  by  the  relative  flow  over  these 
convective  obstacles,  a  process  analogous  to  orographic  grav¬ 
ity  wave  generation  (Clark  et  al.,  1986).  This  mechanism  has 
been  identified  in  observations  of  waves  above  the  convective 
boundary  layer  (Kuettner  et  al.,  1987).  It  has  also  been 
applied  to  time-dependent  clouds  as  a  “transient  mountain” 
mechanism  (Pfister  et  al.,  1993).  It  is  noted  in  passing  that 
when  a  stationary  convective  source  is  considered  there  is  no 
major  difference  between  the  thermal  forcing  and  obstacle 
mechanisms  in  that  the  gravity  waves  are  generated  only 
when  there  is  a  non-zero  background  wind  relative  to  the 
forcing  (Fovell  et  al.,  1992;  Chun  and  Baik,  1998).  When  a 
transient  source  is  considered,  however,  gravity  waves  can  be 
generated  even  for  zero  storm-relative  wind,  and  thus  the  two 
mechanisms  are  distinguishable  in  this  case. 

The  third  mechanism  is  the  so-called  “mechanical  oscilla¬ 
tor”  (Clark  et  al.,  1986;  Fovell  et  al.,  1992;  Lane  et  al.,  2001). 
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Here,  oscillatory  motions  within  a  convective  system  impinge 
on  the  tropopause  and  generate  gravity  waves  in  stable 
regions  above  the  clouds  under  zero  background  wind  relative 
to  the  convective  system.  As  noted  by  Fritts  and  Alexander 
(2003),  this  mechanism  is  very  similar  to  a  diabatic/thermal 
forcing  mechanism  in  which  the  period  of  the  thermal  forcing 
term  is  fixed.  It  is  also  similar  to  an  “oscillatory  obstacle” 
forcing  process  (e.g.,  Prusa  et  al.,  1996). 

Even  though  these  three  mechanisms  are  frequently  dis¬ 
cussed  and  investigated,  a  fundamental  question  exists  as  to 
whether  these  mechanisms  are  really  distinguishable  or  are 
just  different  ways  of  looking  at  the  same  mechanism.  This 
issue  has  not  been  answered  clearly  yet.  Lane  et  al.  (2001) 
calculated  source  terms  of  the  linear  gravity  waves  induced 
by  a  mesoscale  convective  system,  and  showed  that  the  mag¬ 
nitudes  of  non-linear  momentum  and  heat  advection  terms 
(representing  mechanical  oscillator  forcing)  are  much  larger 
than  the  diabatic  heating  term.  From  this  result,  they  conclud¬ 
ed  that  the  mechanical  oscillator  mechanism  dominates  the 
thermal  forcing  mechanism.  However,  recent  2D  numerical 
modelling  by  Song  et  al.  (unpublished  manuscript,  2003) 
reveals  that  the  amplitudes  of  convectively  forced  gravity 
waves  in  the  stratosphere  are  largely  determined  by  the  wave 
propagation  conditions  defined  within  the  horizontal 
wavenumber-frequency  domain  for  given  basic-state  wind 
and  stability  profiles,  rather  than  by  the  magnitude  of  wave 
forcing  terms.  Song  et  al.  (unpublished  manuscript,  2003) 
performed  quasi-linear  simulations  forced  separately  by  non¬ 
linear  advection  and  diabatic  sources  in  the  troposphere,  and 
showed  that  although  the  magnitude  of  the  non-linear  advec¬ 
tion  source  is  much  larger  than  that  of  the  diabatic  source, 
similar  to  the  result  of  Lane  et  al.  (2001),  the  gravity  waves 
radiated  into  the  stratosphere  were  comparable  to  each  other 
in  both  their  amplitude  and  spectral  characteristics.  They 
went  on  to  show  that  vertical  propagation  conditions  restrict¬ 
ed  the  effectiveness  of  some  wave  sources,  such  that  a  large 
portion  of  the  non-linear  advection  forcing  could  not  produce 
gravity  waves  in  the  stratosphere. 

More  detailed  reviews  of  convective  wave  generation  are 
given  in  Section  3  of  Garcia  (2000)  and  in  Sections  3.1.2  and 
7.1  of  Fritts  and  Alexander  (2003).  Detailed  2D  and  3D 
mesoscale  modelling  has  provided  new  insights  into  gravity 
waves  generated  by  realistic  convective  systems  (Hauf  and 
Clark,  1989;  Alexander  et  al.,  1995;  Piani  et  al.,  2000;  Lane 
et  al.,  2001).  From  these  studies  it  has  emerged  that  the  waves 
radiated  from  convection  seem  to  be  interpretable  using  linear 
gravity  wave  theory  (e.g.,  Alexander,  1996;  Pandya  and 
Alexander,  1999;  Lane  et  al.,  2001),  which  greatly  simplifies 
their  parametrization. 

There  have  been  several  attempts  to  parametrize  GWD 
induced  by  subgrid-scale  cumulus  convection  in  large-scale 
models  using  different  theoretical  models  for  determining 
cloud-top  wave  stresses.  Lindzen  (1984)  assumed  that  gravity 
wave  momentum  flux  is  proportional  to  the  square  of  vertical¬ 
ly  integrated  convective  mass  flux,  and  this  parametrization 
was  tested  in  the  global  model  simulations  of  Rind  et  al.  (1988). 


A  later  parametrization  by  Kershaw  (1995)  assumed  that  the 
wave  energy  above  the  convective  region  is  proportional  to  the 
convective  kinetic  energy  in  the  cloud  region.  Bosseut  et  al. 
(1998)  used  a  scheme  in  a  climate  model  where  wave  momen¬ 
tum  fluxes  were  proportional  to  the  precipitation  flux  from  the 
convective  scheme,  with  beneficial  impacts.  Chun  and  Baik 
(1998)  derived  an  analytical  expression  for  the  gravity  wave 
momentum  flux  forced  by  specified  stationary  diabatic  forcing 
representing  the  latent  heating  by  cumulus  convection  for  uni¬ 
form  background  wind  and  buoyancy  frequency.  This  formula¬ 
tion  of  the  cloud-top  wave  stress  has  been  found  to  be 
analogous  to  that  of  surface  mountain  drag  (Pierrehumbert, 
1986)  if  the  non-dimensional  mountain  height  is  replaced  by 
the  “non-linearity  factor”  of  thermally  induced  internal  gravity 
waves  (Lin  and  Chun,  1991).  The  cloud-top  wave  stress  was 
shown  to  be  inversely  proportional  to  the  basic-state  wind 
speed  and  the  cube  of  the  buoyancy  frequency  for  a  given  dia¬ 
batic  heating  rate.  The  formulation  by  Chun  and  Baik  (1998) 
has  been  implemented  in  an  atmospheric  GCM  (Chun  et  al., 
2001)  and  also  in  a  global  spectral  forecast  model  (Kim  and 
Hogan,  unpublished  manuscript,  2003).  This  parametrization 
has  also  been  updated  by  including  basic-state  wind  shear  in  the 
convective  region  and  stability  differences  between  the  con¬ 
vective  region  and  the  regions  above  it  (Chun  and  Baik,  2002). 

In  addition  to  its  impact  in  the  tropics,  the  parametrization 
of  convective  GWD  may  offer  the  additional  benefit  of  alle¬ 
viating  model  climate  biases  in  the  southern  hemisphere 
(Chun  et  al.,  2001)  where  orographic  GWD  has  much  more 
limited  impact  due  to  fewer  mountains.  Furthermore,  in  the 
extratropical  summer  middle  atmosphere,  orographic  gravity 
waves  cannot  penetrate  far  vertically  due  to  critical  levels  (see 
Fig.  1),  and  thus  the  reversal  of  the  summer  stratopause  jets 
must  be  driven  by  non-orographic  GWD  (e.g.,  Jackson, 
1993).  This  issue  is  expanded  upon  in  the  next  section. 

4  Treatment  of  the  middle  atmosphere  in  gravity  wave 
drag  parametrization 

Most  earlier  and  many  current  NWP  models  do  not  include  a 
middle  atmosphere  (heights  -20-90  km).  Possible  reasons  for 
this  are  numerous;  it  adds  substantially  to  the  computational 
burden;  the  existence  of  the  middle  atmosphere  can  be  con¬ 
sidered  unnecessary  for  “short-”  and  “medium-range”  fore¬ 
casting;  and  there  may  be  insufficient  data  to  assimilate  in  the 
middle  atmosphere.  An  extension  toward  “long-range”  fore¬ 
casts,  however,  requires  the  inclusion  of  the  middle  atmos¬ 
phere,  which  affects  the  lower  atmosphere  through  the 
“downward  control”  induced  by  upper-level  drag  (Haynes 
et  al.,  1991)  and  other  dynamical  coupling  effects  (e.g., 
Ambaum  and  Hoskins,  2002).  Indeed  it  has  been  argued  that, 
for  periods  greater  than  one  week,  it  is  necessary  to  know  the 
initial  state  of  the  entire  global  atmosphere  from  the  stratos¬ 
phere  to  the  surface,  as  well  as  the  state  of  the  upper  layers  of 
the  oceans  (Holton,  1992,  based  on  Smagorinsky,  1967).  An 
emerging  body  of  work  is  the  verification  of  the  importance 
of  the  middle  atmosphere  for  long-range  forecasting  of  win¬ 
ter  weather  in  the  northern  hemisphere  (Baldwin  and 
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Dunkerton,  2001;  Thompson  et  al.,  2002).  Global  NWP  mod¬ 
els  are  also  including  a  middle  atmosphere  to  improve  the  for¬ 
ward  modelling  of  satellite  radiances  in  the  operational  data 
assimilation  procedures  that  provide  initial  atmospheric  states 
for  the  forecast  model. 

Thus,  there  is  a  trend  to  extend  the  vertical  domains  of  3D 
large-scale  models  (Pawson  et  al.,  2000).  As  these  upper 
boundaries  of  NWP  and  climate  models  are  extended  into  the 
middle  and  upper  atmosphere,  GWD  parametrization 
becomes  more  important  and  more  challenging.  Since  the 
winds  are  generally  stronger  at  high  altitudes,  the  model  time 
steps  must  be  decreased  for  numerical  stability  thereby  limit¬ 
ing  the  increase  of  spatial  resolution  to  maintain  computa¬ 
tional  efficiency.  The  stratosphere  and  mesosphere  are 
impacted  much  more  profoundly  by  GWD  than  the  tropos¬ 
phere  (e.g.,  Hamilton,  1996;  Shepherd,  2000).  Since  unsatu¬ 
rated  gravity  wave  momentum  flux  densities,  p  u'w',  are 
conserved  in  the  absence  of  dissipation  (see  Appendix  Adi 
for  background),  the  exponential  decrease  of  atmospheric 
densities  p  with  height  means  that  gravity  wave  velocity  fluc¬ 
tuations  [e.g.,  u(x,z,t)]  tend  to  grow  exponentially  in  ampli¬ 
tude  with  height  (Lindzen,  1981).  Thus,  waves  with  miniscule 
amplitudes  in  the  lower  atmosphere  can  attain  large  breaking 
amplitudes  at  upper  levels.  Basic  climatological  wind  patterns 
throughout  the  global  mesosphere  (z~50-90  km)  and  the 
entire  equatorial  middle  atmosphere  (z~  I 5-90  km)  cannot  be 
reproduced  without  parametrizing  GWD  at  these  altitudes 
(McLandress,  1998;  Holton  and  Alexander,  2000). 

In  this  section,  a  general  overview  of  the  current  status  of 
GWD  parametrization  for  the  middle  atmosphere  is  provided. 
More  exhaustive  reviews  can  be  found  elsewhere:  an  intro¬ 
duction  to  some  of  the  single-wave  and  spectral  parametriza- 
tions  for  the  middle  atmosphere  is  given  in  Section  5  of 
McLandress  (1998);  a  collection  of  articles  on  gravity  wave 
parametrization  issues  and  their  performance  in  middle 
atmospheric  models  is  provided  in  the  volume  edited  by 
Hamilton  (1997a);  Fritts  and  Alexander  (2003)  provide  a 
thorough  overview  of  single-wave  and  spectral  gravity  wave 
theories  and  their  transition  to  parametrization  schemes.  A 
historical  perspective  is  provided  by  Hamilton  (1999). 

a  Critical  Level  Filtering  of  Gravity  Waves 
Gravity  waves  generally  dissipate  totally  and  are  absorbed 
into  the  flow  as  they  approach  their  critical  level  zc  where 
their  intrinsic  horizontal  phase  speed  and  vertical  wavelength 
vanish  (as  discussed  in  Appendices  Ac3  and  Ad4).  This  “crit¬ 
ical  level  filtering”  of  gravity  waves  is  important  in  middle 
atmospheric  GWD  parametrization. 

In  a  2D  configuration,  a  hydrostatic  irrotational  gravity 
wave  of  ground-based  horizontal  phase  speed  c  has  a  critical 
level  zc  where  the  condition  c  =  U(zc)  is  met,  where  U(z)  is 
the  background  wind  profile.  Since  wave  phase  speeds  c  are 
constant  in  GWD  parametrizations,  one  can  easily  inspect 
wind  profiles  U{z)  to  see  which  waves  will  encounter  critical 
levels  at  what  heights.  Figure  1  does  this  graphically  for  a  typ¬ 
ical  mid-latitude  location  during  winter  and  summer. 


During  northern  winter  (Fig.  la),  stationary  orographic 
gravity  waves  (c  =  0)  typically  do  not  encounter  critical  lev¬ 
els  and  thus  can  propagate  through  the  middle  atmosphere.  If 
we  assume  that  convection,  on  average,  tends  to  generate 
waves  with  large  eastward  and  westward  phase  speeds 
(Section  3b),  then  the  waves  with  westward  phase  speeds  (c  < 
0)  also  reach  the  mesosphere.  However,  the  eastward  phase 
speed  waves  (c  >  0)  encounter  critical  levels  due  to  the  strong 
eastward  stratospheric  jet.  Similarly,  any  gravity  waves 
launched  from  the  tropospheric  jet  stream  (see  Section  7h), 
which  here  are  assumed  to  have  phase  speeds  similar  to  those 
of  local  jet  stream  wind  speeds,  also  are  absorbed  at  critical 
levels  in  the  stratosphere.  This  leaves  a  preponderance  of  sta¬ 
tionary  and  westward  phase  speed  waves  in  the  mesosphere, 
where  c  -  U(z)  <  0,  and  thus  wave  breaking  in  this  region  will 
drag  winds  strongly  to  the  west  (as  discussed  in  Appendix 
Ad3).  Note  the  importance  of  critical  level  filtering:  without 
the  removal  of  eastward  phase  speed  waves,  this  westward 
drag  would  be  largely  cancelled  by  a  compensating  eastward 
drag  from  c  >  0  waves.  Instead,  these  waves  dissipate  lower 
down  at  stratospheric  critical  levels,  and  should  accelerate 
winds  there  slightly:  recent  work  suggests  that  this  lower 
stratospheric  non-stationary  GWD  may  be  significant  for 
stratospheric  climate  (Alexander  and  Rosenlof,  1996).  Note 
too  that  in  the  lower  stratosphere,  orographic  intrinsic  phase 
speeds  become  small  and  thus  orographic  gravity  waves  tend 
to  saturate  there  as  well  (see  Appendix  Ad),  producing  some 
additional  westward  orographic  GWD  (Lindzen,  1985; 
Palmer  et  al.,  1986;  McFarlane,  1987). 

In  northern  summer  (Fig.  lb)  the  situation  reverses.  Now  a 
strong  westward  stratospheric  jet  absorbs  mountain  waves  and 
westward  phase  speed  gravity  waves,  whereas  the  jet-stream 
waves  and  eastward  phase-speed  waves  from  convection 
reach  the  mesosphere.  This  immediately  highlights  the  impor¬ 
tance  of  non-stationary  GWD  parametrizations  to  the  summer 
extratropical  middle  atmosphere  (e.g.,  Jackson,  1993).  A 
qualitatively  similar  situation  in  the  equatorial  middle  atmos¬ 
phere  makes  non-stationary  waves  from  deep  tropical  convec¬ 
tion  important  all  year  around  (see  also  Sections  3b  and  6a). 

b  Single-Wave  Gravity  Wave  Drag  Parametrization 
Rayleigh  friction  was  used  in  the  first  generation  of  global  3D 
models  with  middle  atmospheres  (Section  2a),  followed  by 
single-wave  Lindzen-type  GWD  parametrization  schemes 
(Section  3a).  The  latter  parametrizations  are  natural  vertical 
extensions  of  the  first  generation  lower  atmospheric  schemes 
and  can  usually  provide  the  orographic  component  of  the 
middle  atmosphere  GWD  in  models  (see  e.g.,  McFarlane  et 
al.,  1997;  Manzini  and  McFarlane,  1998;  Pawson  et  al.,  1998; 
Scaife  et  al.,  2000;  Jackson  et  al.,  2001). 

In  generalized  Lindzen  schemes,  a  collection  of  individual 
waves  with  assigned  phase  speeds  c  are  propagated  vertically 
from  below  under  a  hydrostatic  saturation  model  (Appendix 
Ad2),  assuming  no  interactions  among  the  waves.  These 
schemes  must  be  carefully  tuned  by  varying  wave  parameters 
such  as  amplitudes,  phase  speeds  and  launch  height  at  the 
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source  until  acceptable  results  are  obtained.  Lindzen  schemes 
tend  to  produce  large  rapid  onset  of  drag  with  height,  which 
may  degrade  some  model  simulations  (see  e.g.,  Norton  and 
Thuburn,  1999)  and  may  necessitate  further  tuning  of  para¬ 
meters  or  inclusion  of  some  kind  of  smoothing  (e.g.,  Holton, 
1982;  McLandress  and  McFarlane,  1993;  Hamilton,  1997b). 
A  common  solution  is  to  invoke  an  “intermittency”  factor 
and/or  varying  saturation  thresholds  that  help  scale  down  and 
vertically  redistribute  the  drag,  while  the  factor  and  thresh¬ 
olds  can  also  be  tuned  (Holton,  1982;  Hamilton,  1997b; 
Norton  and  Thuburn,  1999;  Alexander  and  Dunkerton,  1999). 
Despite  the  problems  of  sensitivity  and  arbitrary  tuning,  these 
schemes  can  reproduce  the  basic  circulation  features  of  the 
middle  atmosphere  and,  until  recently,  were  the  standard  way 
of  parametrizing  GWD  at  these  altitudes  in  3D  models  (see 
e.g..  Rind  et  al.,  1988;  Jackson  and  Gray,  1994;  Hamilton, 
1997b;  Norton  and  Thuburn,  1999). 

c  Spectral  Gravity  Wave  Drag  Parametrization 
As  waves  propagate  higher  into  the  atmosphere,  their  wave 
amplitudes  grow  and  propagation  paths  start  to  intersect 
(Fig.  10).  Increased  amplitudes  mean  greater  potential  for 
both  instabilities  (saturation)  and  interactions  among  the 
waves.  Consistent  with  this  idea,  improved  observations  dur¬ 
ing  the  1980s  revealed  broad  frequency  and  wavenumber 
power  spectral  densities  of  gravity  wave-induced  wind  and 
temperature  perturbations  in  the  middle  atmosphere  (see 
reviews  by  Fritts,  1989;  Gardner,  1996;  Section  4  of  Fritts  and 
Alexander,  2003).  It  was  found  that  these  spectra  had  quite 
reproducible  shapes  (VanZandt,  1982;  Smith  et  al.,  1987).  In 
the  case  of  the  vertical  wavenumber  spectrum,  spectral  densi¬ 
ties  at  high  wavenumbers  were  almost  constant  with  height 
and  location,  again  suggestive  of  some  kind  of  saturation 
mechanism  acting  on  the  wave  field.  This  spurred  develop¬ 
ment  of  a  range  of  competing  theories  to  explain  how  these 
intriguingly  reproducible  spectra  are  formed  and  maintained. 
Almost  all  invoke  some  kind  of  “spectral  saturation”  model 
(see  Appendix  Ad2).  Recent  reviews  of  this  area  are  available 
in  the  literature  (Gardner,  1996;  Section  6.3  of  Fritts  and 
Alexander,  2003). 

The  recent  profusion  of  new  spectral  GWD  parametriza- 
tions  has  been  motivated  in  part  by  this  ongoing  debate.  The 
uncertainties  surrounding  it  have  led  to  two  separate 
approaches  to  formulating  spectral  GWD  schemes.  One  class 
of  spectral  parametrizations  retains  close  connections  to 
observed  spectral  shapes,  but  tries  to  simplify  the  spectral  sat¬ 
uration  issue  by  making  basic  assumptions  about  it  when  for¬ 
mulating  drag  algorithms  (Fritts  and  VanZandt,  1993;  Fritts 
and  Lu,  1993a,  1993b;  Warner  and  McIntyre,  1996,  1999, 
2001).  The  approach  is  to  build  single-wave  propagation 
physics  into  a  spectrum  of  component  waves  constrained  in 
amplitude  by  observed  wave  spectra.  These  schemes  then 
generate  GWD  by  saturating  individual  waves  that  reach 
amplitudes  in  excess  of  limits  imposed  by  observed  spectral 
shapes  thereby  preserving  both  basic  single-wave  physics  and 
observed  spectral  distributions  of  wave  energy. 


The  second  class  of  parametrizations  takes  the  opposite 
view:  they  stress  their  own  specific  spectral  theory  of  wave 
saturation  processes  in  generating  a  GWD  parametrization 
(Hines,  1997a,  1997b;  Medvedev  and  Klaassen,  1995,  2000). 
This  approach  tries  to  invoke  wave-wave  interactions  to  pro¬ 
vide  a  more  physical  basis  for  removing  arbitrary  intermitten¬ 
cy  factors  and  artificial  smoothing  of  drag  profiles,  etc.  The 
different  non-linear  wave  physics  in  these  parametrizations 
spawns  different  nomenclature:  “critical-level  filtering”  and 
“wave  breaking”  are  replaced  by  model-dependent  concepts 
such  as  “wave  obliteration”,  “non-linear  wave-diffusive 
damping”,  and  “Doppler  spreading”. 

Another  “spectral”  GWD  parametrization  has  recently 
been  added  to  the  mix  (Alexander  and  Dunkerton,  1999). 
While  advertised  as  spectral,  it  really  differs  from  the  afore¬ 
mentioned  spectral  schemes  in  that  it  makes  no  attempt  to  uti¬ 
lize  or  model  the  observed  shapes  of  gravity  wave  power 
spectral  densities.  In  fact,  it  might  better  be  described  as  a 
“continuous  single-wave”  parametrization  that  improves 
upon  the  “discrete  single-wave”  formulation  of  the  Lindzen- 
like  schemes  discussed  in  Section  4a.  The  latter  schemes 
choose  a  collection  of  discrete  phase  speeds  c  and  propagate 
them  successively  through  the  atmosphere  (Jackson,  1993; 
Hamilton,  1997b;  Norton  and  Thuburn,  1999).  The 
Alexander-Dunkerton  parametrization  uses  a  Lindzen-like 
approach  to  model  wave  propagation  and  saturation, 
but  extends  it  by  incorporating  continuous  phase-speed  distri¬ 
bution  functions,  modifying  an  approach  originally  used  for 
the  QBO  by  Lindzen  and  Holton  (1968).  This  leads  to  differ¬ 
ences  in  the  GWD  calculation:  waves  dissipate  totally  at  their 
saturation  altitude,  rather  than  partially  as  in  the  Lindzen 
(1981)  model. 

These  spectral  parametrizations  also  offer  some  attractive 
practical  advantages  in  models.  The  GWD  profiles  they  pro¬ 
duce  tend  to  be  smoother  than  those  typically  generated  by 
Lindzen-like  models  (Alexander  and  Dunkerton,  1999), 
which  reduces  the  possibility  of  model  instabilities.  The 
observational  constraint  of  spectral  shapes  should  also  mean 
that  they  demand  less  arbitrary  tuning  than  Lindzen-like 
parametrizations.  Nonetheless,  they  all  require  source  specifi¬ 
cations  as  a  lower  boundary  condition  which,  given  the  cur¬ 
rent  uncertainties  about  non-stationary  gravity  wave  sources, 
in  practice  implies  arbitrary  parameter  choices  and/or  exten¬ 
sive  tuning. 

Despite  this  range  of  spectral  parametrization  approaches, 
they  all  respond  in  the  same  general  manner  as  single  waves 
in  Fig.  1:  a  spectrum  of  wave  momentum  flux  density  is  dis¬ 
tributed  at  source  levels  among  a  continuous  distribution  of 
wave  phase  speeds  c,  and  is  propagated  vertically  into  the 
atmosphere.  Certain  portions  are  filtered  out  by  the  back¬ 
ground  winds,  and  surviving  portions  of  the  spectrum  saturate 
at  various  heights  to  generate  a  directional  force  that  either 
accelerates  or  decelerates  the  background  flow.  Nonetheless, 
important  differences  emerge  among  the  GWD  profiles  pre¬ 
dicted  by  various  schemes  (e.g.,  Lawrence,  1997a; 
McLandress,  1997;  Osprey,  2001;  Charron  et  al.,  2002). 
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Given  current  uncertainties  in  non- stationary  wave  sources, 
to  date  spectral  GWD  parametrizations  have  been  used  in  3D 
models  mostly  to  parametrize  global  middle  atmospheric 
GWD  from  indistinct  non-orographic  sources  (e.g.,  Mengel  et 
al.,  1995;  Manzini  and  McFarlane,  1998;  Scaife  et  al.,  2000; 
Scinocca,  2002).  They  appear  to  be  promising  tools:  they  var¬ 
iously  generate  equatorial  QBOs  and  SAOs  (Mengel  et  al., 
1995;  Scaife  et  al.,  2000;  Medvedev  and  Klaassen,  2001)  and 
realistic  extratropical  wind,  temperature  and  eddy  diffusion 
climatologies  (McFarlane  et  al.,  1997;  Manzini  and 
McFarlane,  1998;  Medvedev  et  al.,  1998;  Akmaev,  2001; 
Scinocca,  2002). 

5  Evaluation  of  gravity  wave  drag  parametrization 

a  Online  Evaluation 

Any  GWD  scheme  needs  to  be  tested  within  a  climate  or 
NWP  model  to  assess  its  worth.  Accurate  objective  evalua¬ 
tions  of  a  given  scheme  in  a  large-scale  model  are,  however, 
nontrivial. 

The  first  problem  is  that  it  is  impossible  to  find  a  “perfect” 
model  in  which  only  the  subgrid-scale  GWD  processes  are 
not  represented.  Therefore,  it  is  important  that  the  simulated 
atmosphere  without  any  GWD  parametrization  is  well  under¬ 
stood  in  its  own  right.  When  a  GWD  scheme  is  added,  it  may 
accentuate  certain  model  errors  or  act  to  mask  errors  in  other 
parametrizations  or  ad-hoc  drag  terms  (e.g.,  Klinker  and 
Sardeshmukh,  1992;  Milton  and  Wilson,  1996;  Kim  et  al., 
1998).  For  instance,  Rayleigh  friction,  Newtonian  cooling 
and  other  arbitrary  diffusion  or  dissipation  processes  are  often 
used  in  models  to  tune  and/or  stabilize  the  model.  While  all 
this  interacting  damping  can  seem  beneficial  in  terms  of  the 
simulated  monthly  zonal-mean  fields,  it  can  often  suppress 
realistic  dynamical  variability  within  the  model  (Shepherd 
et  al.,  1996;  Lawrence,  1997b;  Pawson  et  al.,  1998).  This 
problem  is  further  exacerbated  by  the  fact  that  most,  if  not  all, 
GWD  schemes  contain  parameters  that  are  not  well  con¬ 
strained  and  must,  in  practice,  be  varied  or  “tuned”  until  the 
simulated  climate  looks  reasonable.  Such  tuning  makes  it  dif¬ 
ficult  to  assess  any  given  GWD  scheme  objectively.  In  some 
cases,  tuning  GWD  is  erroneously  regarded  as  a  final  means 
to  fix  all  the  remaining  biases  in  the  simulated  atmosphere, 
and  such  an  approach  can  lead  to  problems.  For  example,  if 
the  model’s  ozone  heating  parametrization  is  misrepresented 
in  such  a  way  as  to  induce  an  excessively  strong  stratospher¬ 
ic  jet,  tuning  a  GWD  scheme  can  mask  the  problem,  an  effect 
that  clearly  does  not  improve  the  model  (Kim  et  al.,  1998). 

Even  neglecting  such  problems,  there  are  other  issues  to 
consider  when  assessing  model  simulations  performed  with, 
or  without,  GWD.  Observational  assessments  of  GWD 
parametrization  schemes  can  be  performed;  e.g.,  by  a  network 
of  medium  frequency  (MF)  radar  observations  comparing 
measured  winds  and  waves  with  those  generated  by  a  middle 
atmospheric  model  with  parametrized  GWD  (Manson  et  al., 
2002).  With  climate  models,  however,  simulations  are  usual¬ 
ly  performed  for  months  or  years,  and  monthly  and  zonal- 
mean  averages  are  composed  from  model  output  and 


compared  with  observed  climatology  as  a  basic  performance 
indicator  for  a  GWD  scheme.  NWP  models  can  also  be 
assessed  in  “climate  mode”  without  invoking  periodic  data 
assimilation  cycles:  however,  short-  and  medium-range  fore¬ 
cast  integrations  are  usually  performed  using  various  initial¬ 
ization  and  data  assimilation  procedures  which  makes  the 
GWD  scheme’s  net  impact  harder  to  evaluate.  Assessments 
based  on  “skill  scores”  of  standard  metrics  such  as  the  root- 
mean-square  (RMS)  error  and  anomaly  correlation  of  the 
geopotential  heights  can  be  useful.  In  this  case  the  perfor¬ 
mance  of  a  GWD  scheme  is  highly  dependent  on  its  harmony 
with  the  initialization  and  assimilation  schemes  as  well  as  its 
physical  adequacy.  Another  useful  approach  has  been  to  diag¬ 
nose  the  large-scale  momentum  balance  of  the  simulated 
atmosphere  at  various  time  steps  by  outputting  all  the  mod¬ 
elled  and  parametrized  terms  in  the  momentum  balance  equa¬ 
tion  (Bell  et  al.,  1991;  Klinker  and  Sardeshmukh,  1992; 
Milton  and  Wilson,  1996).  “Single  column”  versions  of  the 
model,  which  are  often  used  to  test  convection  and  radiation 
parametrizations  (e.g.,  Wu  and  Moncrieff,  2001),  can  also  be 
used  to  test  GWD  column  parametrizations. 

b  Offline  Evaluation 
1  AGAINST  OBSERVATIONS 

In  addition  to  testing  and/or  tuning  GWD  schemes  within 
large-scale  models,  it  is  useful  to  test  the  schemes  “offline” 
against  actual  observations  of  waves  and  GWD.  Due  to  the 
relative  ease  of  making  low-level  atmospheric  measurements 
over  land,  orographic  GWD  schemes  are  amenable  to  these 
observational  assessments.  Near-surface  mountain  (pressure) 
drag  predictions  (see  Appendix  Aa)  have  been  evaluated  using 
data  from  microbarograph  chains  (e.g..  Smith,  1978;  Vosper 
and  Mobbs,  1997;  Miranda  et  al.,  1999).  More  rigorous  tests 
of  GWD  predictions  have  been  made  using  data  from  multi¬ 
instrument  field  campaigns  over  major  mountain  ranges  (e.g., 
Bougeault  et  al.,  1993,  2001;  Beau  and  Bougeault,  1998). 

Since  convective  GWD  parametrizations  are  more  recent, 
currently  observations  help  most  in  developing  rather  than 
evaluating  these  schemes.  Clear  correlations  between  convec¬ 
tive  activity  and  enhanced  gravity  wave  energy  have  been 
observed,  particularly  in  the  tropics  (e.g.,  Tsuda  et  al.,  1994, 
2000;  Karoly  et  al.,  1996;  Alexander  et  al.,  2000;  McLandress 
et  al.,  2000;  Preusse  et  al.,  2001).  Wave  momentum  fluxes 
and  GWD  have  been  estimated  from  observations  of  waves 
near  convection  in  various  ways  (e.g.,  Sato,  1993;  Vincent 
and  Alexander,  2000;  Alexander  et  al.,  2000). 

While  spectral  models  of  gravity  wave  variances  were  orig¬ 
inally  motivated  (and  constrained)  by  observations,  their 
GWD  predictions  are  harder  to  evaluate  observationally. 
Eckermann  (1995)  has  shown  that  observed  gravity  wave 
spectra,  rather  than  being  “universal”  as  was  originally 
assumed  (VanZandt,  1982),  vary  in  response  to  background 
winds  in  a  manner  consistent  with  some  spectral  GWD 
parametrizations.  Other  observed  spectral  sensitivities  (e.g., 
Nastrom  et  al.,  1997;  Whiteway,  1999)  seem  to  be  at  odds 
with  standard  model  predictions. 
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2  AGAINST  MESOSCALE  MODEL  SIMULATIONS 
While  measurements  should  ideally  be  used  for  validation, 
this  can  prove  difficult  in  practice.  Although  wave  perturba¬ 
tion  profiles  can  be  measured  fairly  easily,  gravity  wave 
momentum  flux,  p  u'w ' ,  the  key  wave  quantity  for  GWD  cal¬ 
culations  (see  Appendices  Aa  and  Ad)  is  notoriously  difficult 
to  measure  directly  (e.g.,  Worthington  and  Thomas,  1996). 
Further,  when  measured  it  proves  to  be  an  intrinsically 
“noisy''  quantity  (Kudeki  and  Franke,  1998),  and  thus  its  ver¬ 
tical  divergence  —  the  GWD  —  can  be  still  noisier  and  thus 
uncertain  (e.g.,  Chang  et  al.,  1997). 

Since  high-resolution  mesoscale  models  can  simulate 
regional  wave  breaking  episodes  explicitly,  GWD  can  be  cal¬ 
culated  from  a  mesoscale  model's  output  and  compared  with 
the  prediction  of  a  simpler  GWD  scheme.  Some  of  the  earli¬ 
est  GWD  parametrizations  were  compared  with  available 
analytic  solutions  based  mainly  on  2D  linear  single  stationary 
waves  over  an  isolated  idealized  mountain  (e.g.,  Wurtele  et 
al.,  1987),  and  then  against  progressively  more  complex 
mesoscale  model  simulations  of  orographic  flows  (Kim  and 
Arakawa,  1995;  Broad,  1996;  Flereil  and  Stein,  1999)  and 
convective  systems  (e.g.,  Kershaw,  1995). 

Operational  3D  mesoscale  models  can  “hindcast”  particu¬ 
lar  wave  events  in  great  detail  (e.g..  Powers,  1997;  Doyle  and 
Shapiro,  2000),  and  so  can  provide  a  useful  link  between 
imperfect  data  and  GWD  schemes.  For  validation,  a  detailed 
3D  mesoscale  model  simulation  of  a  measured  wave  event 
can  be  performed,  and  its  results  can  be  compared  with  avail¬ 
able  data.  If  those  comparisons  are  acceptable,  this  “validat¬ 
ed”  mesoscale  model  output  can  be  used  to  compute  the 
GWD,  which  can  then  be  compared  to  that  predicted  by  the 
GWD  parametrization  (Broad,  1996,  2000).  Recent  high-res- 
olution  global  models  have  resolved  gravity  wave  motions 
whose  spectral  shapes  agree  with  both  observations  and  forms 
implemented  in  spectral  GWD  schemes  (Koshyk  et  al.,  1999). 
The  next  section  discusses  resolved  gravity  waves  in  large- 
scale  models  more  fully. 

6  Resolved  gravity  waves  versus  parametrized  gravity 
wave  drag 

Since  GWD  parametrizations  are  designed  to  act  primarily  on 
synoptic-scale  flows,  they  derive  and  use  equations  that 
assume  that  GWD  acts  on  a  large-scale  quasi-stationary 
atmosphere  largely  devoid  of  any  resolved  gravity  waves.  It 
is  well  known  that,  while  unbalanced  initialization  fields  gen¬ 
erate  spurious  gravity  waves  in  global  models,  balanced  ini¬ 
tialization  procedures  largely  eliminate  these  problems.  Thus, 
these  simplifying  premises  are  appropriate  when  using  these 
parametrizations  in  low-to-medium  resolution  global  models. 

Currently,  however,  global  models  are  being  run  at  suffi¬ 
ciently  fine  spatial  resolutions  that  “real”  gravity  waves  are 
being  generated  explicitly  by  sources  within  the  models. 
Spectral  analysis  of  high-resolution  model  output  in  the  upper 
troposphere  and  lower  stratosphere  based  on  spherical  har¬ 
monic  expansions  reveals  a  transition  in  behaviour  at  total 


spherical  harmonic  wavenumbers  n  ~  80,  from  rotationally 
dominated  flow  at  scales  n  <  80  to  shallower  spectra  at n  >  80 
where  divergent  gravity  wave  motions  are  as  energetic  as,  or 
more  energetic  than,  vortical  motion  (Koshyk  et  al.,  1999; 
Koshyk  and  Hamilton,  2001).  This  breakpoint  wavenumber 
decreases  with  altitude  as  these  gravity  waves  increase  in 
amplitude  with  height.  Since  current  state-of-the-art  spectral 
NWP  models  run  anywhere  in  the  range  -T200-T500  (i.e., 
triangular  spectral  truncation  at  a  total  wavenumber 
-200-500,  implying  horizontal  resolutions  of  Ax  ~0.25°-0.6°, 
i.e.,  -40-100  km),  it  is  clear  that  these  models  are  generating 
significant  amounts  of  gravity  wave  activity  and  that  these 
waves  likely  dominate  small-scale  model  variability  at  upper 
levels.  Thus  it  is  very  important  to  study  these  waves,  partic¬ 
ularly  the  way  in  which  they  contribute  to  GWD  and  interact 
with  any  subgrid-scale  GWD  parametrizations. 

a  Tropical  Gravity  Waves  Generated  by  Convection 
Since  /  — »  0  in  the  tropics,  long-period  gravity  waves  with 
long  zonal  wavelengths  can  exist  there  (see  Appendix  Ac),  as 
well  as  other  tropical  modes  such  as  Kelvin  waves  and  mixed 
Rossby-gravity  waves.  Deep  tropical  convection  generates  a 
rich  spectrum  of  such  waves  in  global  models  run  at  quite 
moderate  resolutions.  For  example,  Ricciardulli  and  Garcia 
(2000)  compared  spectra  of  the  vertical  component  of  the 
Eliasson-Palm  (E-P)  flux  (essentially  p  u'w '  for  gravity 
waves;  Appendix  Adi)  of  tropical  waves  generated  by  vari¬ 
ous  convective  adjustment  parametrizations  in  the 
Community  Climate  Model  (CCM3)  model  (T42L18, 
Av-2. 8°)  with  those  calculated  from  observed  proxies  for 
deep  convective  activity  derived  from  high-resolution  global 
cloud  imagery  (GCI).  Nissen  et  al.  (2000)  conducted  similar 
experiments  with  a  T21  GCM  (Av~5.6°),  while  Amodei  et  al. 
(2001)  compared  tropical  wave  properties  from  a  variety  of 
large-scale  models  (Ax~2.5°-5.5°).  Mesoscale  activity  in 
models  appears  to  be  enhanced  generally  in  the  tropics  rela¬ 
tive  to  the  extratropics  (Koshyk  and  Hamilton,  2001)  in 
agreement  with  observations  (e.g.,  Alexander  et  al.,  2002). 

Ricciardulli  and  Garcia  (2000)  found  large  discrepancies 
between  wave  fields  inferred  from  observed  GCI  convection 
and  those  generated  by  some  convective  schemes  in  models; 
the  latter  significantly  underestimated  wave  momentum  flux¬ 
es,  particularly  gravity  wave  fluxes,  as  shown  in  Fig.  7.  They 
noted  that  convection  schemes  in  models  have  been  devel¬ 
oped  and  tested  mostly  for  their  ability  to  generate  realistic 
time-mean  distributions  of  tropical  convection,  and  suggested 
attention  also  be  paid  to  their  ability  to  reproduce  temporal 
variability  in  convection.  This  assumes  that  the  algorithms 
used  by  Ricciardulli  and  Garcia  (2000)  to  convert  observed 
GCI  data  into  latent  heating  rates  yield  accurate  radiated  wave 
fields.  In  assessing  and  reviewing  various  GCM  studies, 
Amodei  et  al.  (2001)  noted  that  simple  moist  convective 
adjustment  seemed  to  generate  greater  and  more  realistic 
tropical  wave  fields  than  some  of  the  more  sophisticated 
schemes  currently  used  in  models.  These  findings,  while 
interesting,  are  also  preliminary  and  more  work  is  needed  to 
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Fig.  7  Zonal  wavenumber  spectra  PEP(k)  of  vertical  component  of  EP-flux  for  all  equatorial  gravity  waves  with  periods  >  6  hours  generated  by  deep  tropical 
convection  between  15°S-15°N,  derived  by  adding  relevant  component  curves  in  Figs  13b  and  14b  of  Ricciardulli  and  Garcia  (2000)  and  plotting  in 
“energy  content”  or  “area  preserving”  form  kPEp(k).  The  black  curve  is  derived  from  convective  heating  specified  by  3-hourly  global  cloud  imagery 
(GCI)  data.  The  grey  curve  is  output  from  the  CCM3  global  model  using  the  convective  heating  parametrization  of  Hack  (1994).  Curves  are  averaged 
and  smoothed  in  wavenumber  space:  see  Ricciardulli  and  Garcia  (2000)  for  full  details. 


understand  and  evaluate  tropical  convection  and  waves  in 
models,  particularly  with  regard  to  their  roles  in  the  simula¬ 
tion  of  tropical  circulation  features  like  the  QBO  and  SAO. 

b  Sensitivity  to  Vertical  and  Horizontal  Resolutions 
Increasing  the  spatial  resolution  of  a  model  should  obviously 
help  it  resolve  more  gravity  waves,  their  dissipation  and 
GWD.  However,  Lindzen  and  Fox-Rabinovitz  (1989)  argued 
for  consistency  in  choices  of  the  ratio  of  horizontal  to  vertical 
resolution  in  models,  and  claimed  that  excess  horizontal  res¬ 
olution  due  to  inadequately  low  vertical  resolution  might 
increase  model  “noise”.  The  precise  dependencies  of  models 
to  changes  in  vertical  and  horizontal  resolution  are  not  clear 
(e.g..  Lander  and  Hoskins,  1997;  Davies  and  Brown,  2001) 
and  so  modellers  have  explored  these  sensitivities  experimen¬ 
tally  (e.g.,  Hamilton  et  al.,  1999). 

Resolved  tropical  wave  fields  prove  most  sensitive  to 
increases  in  model  vertical  resolution,  Az  (e.g.,  Nissen  et  al., 
2000).  The  consequent  improvements  in  explicitly  simulated 
tropical  wave  fields  and  wave  drag  are  now  cited  as  reasons 
for  the  elusiveness,  until  recently,  of  simulating  a  QBO  in  a 
GCM  without  parametrized  GWD  (see  Section  3b).  While 


models  show  some  sensitivity  to  increased  horizontal  resolu¬ 
tion,  Ax,  this  is  more  of  a  secondary  effect  and  may  be  an  indi¬ 
rect  manifestation  of  other  more  pertinent  influences,  such  as 
horizontal  diffusivity  (e.g.,  Hamilton  et  al.,  2001). 

Almost  the  opposite  dependencies  arise  in  the  extratropics. 
Since  ^=0,  extratropical  gravity  waves  typically  have  shorter 
periods  and,  according  to  the  dispersion  relation  (e.g..  Gill, 
1982;  see  Eq.  (2)  in  Appendix  A),  shorter  horizontal  wave¬ 
lengths  than  in  the  tropics.  Thus  models  must  be  run  at  high¬ 
er  horizontal  resolutions  before  explicitly-resolved  gravity 
waves  appear  prominently  in  the  extratropics.  As  shown  in 
Fig.  8,  improved  horizontal  resolution  yields  larger  and  more 
variable  extratropical  gravity  wave  momentum  fluxes.  At 
Av~l°,  a  fairly  rich  spectrum  of  extratropical  gravity  waves 
emerges  in  most  models  (e.g.,  Hayashi  et  al.,  1989;  Sato  et  al., 
1999).  However,  these  wave  fields  show  much  less  sensitivi¬ 
ty  to  increases  in  vertical  resolution  (e.g.,  O’Sullivan  and 
Dunkeiton,  1995).  Such  dependencies  are  reflected  in  the 
resolved  GWD:  simulations  of  the  extratropical  middle 
atmosphere  without  fully  parametrized  GWD  improve  notice¬ 
ably  when  horizontal  resolution  is  increased  (i.e.,  Ax  is 
reduced),  but  show  much  less  sensitivity  to  increased  vertical 
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Fig.  8  Time  series  of  the  momentum  flux  term  divided  by  pressure,  -  w'coTp,  at  ~60°N  and  0.1  hPa  (~65  km)  for  wavenumbers  k  >  5  in  GFDL  SKYHI  simu¬ 
lations  for  20-28  December.  Here  co'  is  the  vertical  pressure  velocity  perturbation.  The  black  curve  shows  the  results  from  N90  resolution  model  (Ax~l° 
-  1.2°),  the  grey  curve  results  from  N30  resolution  model  (Ax~3°  -  3.6°).  (Adapted  from  Fig.  5  of  Hayashi  et  al.,  1989;  with  permission  from  Birkhauser 
Publishing  Inc.) 


resolution  (Boville,  1991,  1995;  Hamilton  et  al.,  1995,  1999; 
Jones  et  al.,  1997). 

c  An  End  to  GWD  Parametrization? 

Unlike  climate-chemistry  models,  which  still  run  at  fairly 
moderate  resolutions,  some  global  NWP  models  now  equal  or 
even  exceed  the  resolution  of  the  modelling  studies  discussed 
in  the  previous  subsection.  Since  such  models  are  now  simu¬ 
lating  some  GWD,  it  suggests  that  as  resolutions  improve, 
GWD  parametrization  might  soon  become  unnecessary  in 
high-resolution  models.  However,  several  factors  raised  in  the 
following  text  suggest  that  the  spatial  resolutions  necessary  to 
achieve  this  are  still  much  shorter  than  the  highest  resolutions 
attained  to  date  by  global  models,  indicating  that  GWD 
parametrization  will  remain  important  for  both  global  NWP 
and  climate  models  for  some  time  to  come. 

1  RESOLVING  THE  WAVE  OSCILLATIONS 

Gravity  wave  energy  spectra  vary  with  zonal  wavenumber  k 
roughly  as  k~5/i  and  thus  the  bulk  of  the  energy  variance  can 
be  contained  at  larger  resolved  scales  (Hamilton,  1996; 
Koshyk  et  al.,  1999).  On  the  other  hand,  as  pointed  out  by 
Hamilton  (1996),  spectra  of  wave  momentum  flux,  the  most 
relevant  quantity  for  GWD  calculations,  decrease  much  more 
gradually  with  increasing  k.  Figure  9  plots  a  cospectrum  of 
zonal  and  vertical  (pressure)  velocity  perturbations  from  a 
Geophysical  Fluid  Dynamics  Laboratory  (GFDL)  SKYHI 
N90L40  (Ay~1°-1.2°)  simulation  at  -45 °N  during  summer  at 
a  height  of  -0.03  hPa  (Hamilton,  1996).  The  spectrum  is  plot¬ 
ted  in  “energy  content”  form,  so  that  the  area  beneath  it  is  pro¬ 
portional  to  the  total  u  co '  flux  (a  close  proxy  for  momentum 
flux  p  u'w' ),  where  co'  is  the  vertical  pressure  velocity  pertur¬ 
bation.  Beyond  k  ~  30  we  see  that  the  curve  (and  hence  the  flux 
contribution)  increases  gradually  but  monotonically  up  to  the 
smallest  resolvable  scales  in  the  model.  While  these  spectral 


shapes  vary  considerably  (note  that  u  co '  is  signed  and  thus  its 
cospectrum  can  be  positive  or  negative).  Fig.  9  nonetheless 
highlights  the  importance  of  the  smallest  resolved  wave  scales 
to  the  total  wave  momentum  flux.  Observations  and  models 
suggest  that  the  final  break  in  these  spectra,  at  which  the  hori¬ 
zontal  spectrum  of  gravity  wave  oscillations  would  be  fully 
resolved,  occurs  at  scales  0(1-10  km)  (e.g.,  Bacmeister  et  al., 
1996;  Broad,  2000).  This  implies  ~10-100-fold  increases  in 
horizontal  resolution  beyond  even  these  state-of-the-art  simu¬ 
lations  that  already  stretch  current  computing  capabilities. 

2  RESOLVING  THE  SATURATION  PROCESS 
The  previous  arguments  relate  to  resolving  the  wave  oscilla¬ 
tions  themselves.  Wave  saturation  and  GWD,  however, 
involve  “subwavelength-scale”  instabilities.  While  the 
exact  details  are  uncertain  (see  Section  7a  and  Appendix  Ad), 
saturation  generally  involves  formation  of  unstable  turbulent 
layers  that  are  significantly  shorter  fractions  of  the  wave¬ 
lengths  of  the  dissipating  waves  (see,  e.g.,  Lindzen,  1988; 
McIntyre,  1989;  Hines,  1991,  1992).  Thus,  resolving  GWD  in 
models  would  seem  to  require  resolving  these 
unstable  portions  of  a  wave  field,  then  invoking  the  model's 
turbulent  diffusion  parametrization  to  mix  them  and 
thereby  explicitly  saturate  the  waves  (e.g.,  Liu  et  al.,  1999). 
Observations  suggest  these  layers  have  vertical 
widths  <50-200  m  in  the  troposphere  and  lower  stratosphere, 
increasing  perhaps  to  <500-1000  m  in  the  upper  stratosphere 
and  lower  mesosphere  (Sato  and  Woodman,  1982;  Cot 
and  Barat,  1986;  Yamamoto  et  al.,  1987;  Fritts  and  Werne, 
2000),  with  theory  implying  horizontal  widths  -10-100  times 
longer  (e.g.,  Lindzen,  1988).  Such  resolution  in  global 
3D  models  is  well  beyond  current  computing  capabilities: 
further,  it  probably  also  requires  new  generations  of  non¬ 
hydrostatic  global  models  (e.g.,  Daley,  1988;  Smolarkiewicz 
et  al.,  2001). 
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Fig.  9  Zonal  wavenumber  cospectrum  Plia(k)  of  perturbations  of  zonal  velocity  u  and  vertical  pressure  velocity  co'  due  to  gravity  waves  explicitly  resolved  at 
0.03  hPa  (~73  km)  and  ~45°N  in  a  GFDL  SKYHI N90L40  (A.r~l°  -  1.2°)  model  run  for  June/July.  The  result  is  taken  from  Fig.  18  of  Hamilton  (1996), 
with  permission  from  Elsevier  Science  and  replotted  in  "energy  content”  form  kP  (k)  appropriate  in  log-log  space.  The  spectrum  is  an  average  of  many 
individual  spectra  and  is  smoothed  in  wavenumber  space.  The  grey  curve  shows  the  k+l12  power  law. 


Given  this,  it  is  natural  to  ask  how  realistic-looking  GWD 
can  have  been  explicitly  simulated  within  some  current  high- 
resolution  global  middle  atmosphere  models  (see,  e.g.,  Jones 
et  al.,  1997;  Hamilton  et  al.,  1999).  Rather  than  being  saturat¬ 
ed  by  small-scale  turbulent  instabilities  within  the  wave,  it 
seems  likely  that  waves  within  these  models  are  refracted  to 
shorter  wavelengths  and  then  removed  or  “quenched”  by 
large-scale  (i.e.,  wavelength-scale  and  above)  turbulent  and 
numerical  diffusion  schemes  within  the  model  (see,  e.g., 
Fig.  6  of  Marks  and  Eckermann,  1995).  This  has  interesting 
similarities  to  the  total  breakdown  assumption  used  in  some 
GWD  schemes,  which  also  seem  to  produce  reasonable  esti¬ 
mates  of  GWD  (e.g.,  Alexander  and  Dunkerton,  1999).  Thus, 
while  the  precise  saturation  mechanisms  may  be  wrong,  such 
“quenched  GWD”  simulated  in  models  may  be  entirely  ade¬ 
quate  for  current  global  modelling  purposes,  which  might 
then  require  the  comparatively  more  moderate  increases  in 
resolution  implied  in  the  previous  sub-subsection  simply  to 
resolve  the  waves  themselves. 

7  Further  issues  in  gravity  wave  drag  parametrization 

The  discussion  in  the  previous  sections  has  made  clear  that 
GWD  parametrization  for  NWP  and  climate  models  is  an 


important  ongoing  endeavour  that  is  by  no  means  solved.  As 
evident  from  Section  6c,  even  optimistic  projections  of  spatial 
resolution  increases  in  models  do  not  offer  any  imminent 
prospect  of  eliminating  the  need  to  parametrize  GWD. 

Thus,  further  developments  are  essential.  This  section  rais¬ 
es  a  series  of  issues  (many  interrelated)  for  GWD  parame¬ 
trization  that  we  believe  deserve  further  study. 

a  Wave  Saturation 

As  has  been  discussed  (Sections  4b,  4c  and  6c;  Appendix 
Ad),  wave  instabilities  that  lead  to  saturation  and  GWD  are 
still  poorly  understood,  and  so  their  specification  in  GWD 
schemes  remains  fairly  basic. 

While  most  single-wave  schemes  still  use  a  saturation 
model  similar  to  that  of  Lindzen  (1981),  a  number  of  issues 
associated  with  this  have  been  explored  and  extensions  pro¬ 
posed  (see  e.g.,  the  review  of  Dunkerton,  1989).  GWD  and 
vertical  mixing  of  heat  and  momentum  prove  acutely  sensi¬ 
tive  to  uncertainties  in  the  degree  of  wave  supersaturation  and 
turbulence  localization  (Fritts  and  Dunkerton,  1985;  Lindzen, 
1988;  McIntyre,  1989;  Kim  and  Mahrt,  1992;  Liu,  2000).  It 
has  also  been  shown  that  inertia  gravity  waves  break  at  small¬ 
er  amplitudes  through  dynamical  rather  than  convective 
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instabilities  (e.g.,  Lelong  and  Dunkerton,  1998).  Mechanisms 
of  saturation  and  GWD  due  to  the  “resonance  breaking”  of 
trapped  mountain  waves  also  remain  somewhat  uncertain 
(e.g.,  Georgelin  and  Lott,  2001).  The  uncertainties  in  spectral 
saturation  theories  are  even  larger  (Gardner,  1996;  Fritts  and 
Alexander,  2003).  As  discussed  in  Section  4c,  some  spectral 
schemes  assume  that  these  saturation  details  are  relatively 
unimportant  for  GWD  calculations,  while  others  assume  that 
they  are  critical.  All  these  issues  merit  further  study. 

Numerical  studies  suggest  that  much  different  dynamics 
govern  breakdown/saturation  of  gravity  waves  in  3D  than  in 
2D.  Fritts  and  Werne  (2000)  review  some  of  this  recent  work, 
but  note  that  it  is  perhaps  too  early,  at  present,  to  assess  the 
general  impact  for  GWD  and  mixing  (see  also  Fritts  and 
Alexander,  2003).  3D  effects  are  discussed  more  generally  in 
the  following  subsection. 

b  Three  Dimensionality 

Noticeable  differences  in  many  gravity  wave  properties 
emerge  on  progressing  from  2D  to  3D  models.  Isolated  3D 
sources  radiate  characteristic  3D  gravity  wave  patterns,  such 
as  “ship  waves”  from  circular  mountains  (Smith,  1980; 
Sharman  and  Wurtele,  1983;  Shutts,  1998)  and  circular/coni¬ 
cal  gravity  wave  patterns  from  convection  (Piani  et  al.,  2000; 
Lane  et  al.,  2001).  This  implies  directional  GWD  and  critical 
level  effects  that  are  beginning  to  receive  attention  in  param- 
etrizations  (e.g.,  Bacmeister,  1993;  Shutts,  1995;  Broad, 
1995;  Hines,  1997a;  Gregory  et  al.,  1998;  Scinocca  and 
McFarlane,  2000).  Such  3D  waves  radiate  into  progressively 
greater  horizontal  volumes,  an  effect  that  taxes  the  validity  of 
the  column  parametrization  approach  (discussed  in  the  next 
subsection),  particularly  at  upper  levels.  Moreover,  this  geo¬ 
metrical  spreading  of  wave  action  leads  to  conservative  wave 
amplitude  evolution,  saturation  and  GWD  that  can  differ 
markedly  from  standard  2D  theory  (Shutts,  1998;  Vosper  and 
Mobbs,  1998;  Broad,  1999;  Broutman  et  al.,  2001).  3D 
spreading  tends  to  reduce  wave  momentum  fluxes  and  GWD 
at  upper  levels  compared  to  2D  predictions  (e.g.,  Shutts, 
1998).  Both  lateral  spreading  of,  and  reductions  in,  upper- 
level  GWD  seem  to  be  needed  in  many  models  (see 
Section  7d). 

Important  changes  can  also  occur  at  or  near  the  source. 
Modelling  flow  over  3D  mountains  yields  surface  pressure 
drag  and  radiated  wave  momentum  fluxes  that  again  are 
smaller  than  for  2D  ridges  (e.g.,  Clark  and  Farley,  1984; 
Nappo  and  Chimonas,  1992).  Low-level  mountain  wave 
breaking  also  varies  in  form  and  location  on  transitioning 
from  2D  to  3D  obstacles  (e.g..  Smith,  1989;  Bauer  et  al., 

2000) ,  as  does  resonance  breaking  (Miranda  and  Valente, 
1997).  3D  surface  flow  can  be  blocked  upstream,  or  else  pass 
around  the  mountain  where  it  can  generate  vortices  down¬ 
stream,  producing  drag  in  either  case  and  often  reducing  wave 
amplitudes  aloft  (e.g.,  Peng  et  al.,  1995;  Schar  and  Durran, 
1997;  Bauer  et  al.,  2000;  Epifanio  and  Durran,  2001;  Smith, 

2001) .  These  and  other  3D  effects  have  been  partially  imple¬ 
mented  in  some  subgrid-scale  orographic  GWD  schemes 


through  additional  orographic  statistics  (Kim  and  Arakawa, 
1995)  and  are  beginning  to  be  fully  parametrized  in  newer 
GWD  schemes  (e.g.,  Lott  and  Miller,  1997a;  Gregory  et  al., 
1998;  Lott,  1999;  Scinocca  and  McFarlane,  2000)  and  surface 
roughness  schemes  (e.g..  Wood  et  al.,  2001). 

c  Column  Parametrization 

Global  models  prefer  parametrizations  that  operate  in  vertical 
columns  above  each  horizontal  grid  box,  since  they  run  quick¬ 
ly  and  efficiently  on  parallel  architectures.  Thus,  all  existing 
parametrizations  of  GWD,  both  single-wave  (see  Sections  3a, 
3b  and  4b)  and  spectral  (see  Section  4c),  are  formulated  as 
column  algorithms.  Stationary  2D  orographic  gravity  waves 
propagate  almost  vertically,  and  thus  orographic  GWD 
parametrizations  translate  fairly  naturally  into  “column 
parametrizations”.  Conversely,  non-stationary  waves  from 
non-orographic  (and  orographic)  sources  propagate  obliquely 
from  their  source  into  adjacent  horizontal  grid  boxes  (see 
Fig.  10),  as  do  3D  and  non-hydrostatic  orographic  waves.  In 
reality  this  should  yield  a  lateral  spreading  of  GWD  that  has 
been  pointed  out  from  time  to  time  (e.g.,  Kim  and  Arakawa, 
1995;  see  also  previous  subsection).  Nonetheless,  any  param¬ 
etrization  that  attempted  to  model  this  lateral  wave  spreading 
directly  would  probably  add  drastically  to  the  computational 
burden.  Simpler  ways  of  parametrizing  this  may  be  worth 
exploring  since  horizontally  localized  GWD  sometimes 
seems  to  generate  model  errors  (e.g.,  Klinker  and 
Sardeshmukh,  1992). 

d  Low-Level  versus  Upper-Level  Gravity  Wave  Drag 
Parametrization 

As  discussed  in  Section  4b,  single -wave  GWD  schemes  in  the 
middle  atmosphere  tend  to  produce  somewhat  excessive  drag 
rather  suddenly,  and  GWD  profiles  can  be  scaled  down  and/or 
smoothed  out  in  various  ways.  Orographic  GWD  schemes  in 
global  models  often  produce  excessive  “upper-level”  (stratos¬ 
pheric)  GWD  and  insufficient  “low-level”  (lower  tropospher¬ 
ic)  GWD  (Klinker  and  Sardeshmukh,  1992;  Kim,  1996;  Milton 
and  Wilson,  1996;  Kim  et  al.,  1998;  Gregory  et  al.,  1998). 

A  beneficial  impact  has  been  reported  through  an  overall 
ad-hoc  decrease  of  the  ratio  of  upper-level  to  lower-level  oro¬ 
graphic  GWD  in  models  (e.g..  Miller  et  al.,  1989;  Iwasaki 
et  al.,  1989),  or  through  a  selective  and  systematic  decrease  of 
this  ratio  (e.g.,  Kim,  1996;  Milton  and  Wilson,  1996).  Finding 
an  optimal  ratio,  however,  can  be  complicated:  moreover,  the 
ratio  may  not  be  the  only  problem. 

Recent  work  has  focused  on  parametrizing  additional 
sources  of  low-level  GWD  (Section  3a)  and  surface  friction 
drag  (see  Section  7e2),  which  act  to  redistribute  low-level 
drag  in  models  (Kim  and  Hogan,  unpublished  manuscript, 
2003).  Further,  for  orographic  waves,  addition  of  enhanced 
surface  friction  drag  (Section  7e2),  three-dimensionality 
(Section  7b)  and  stagnant  flow  regimes  in  valleys  (e.g..  Smith, 
2001)  all  seem  to  reduce  wave  amplitudes  and  breaking  at 
upper  levels.  Such  work  perhaps  represents  a  new  and  more 
methodical  approach  to  rectifying  this  recurrent  problem. 
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e  Coupling  among  Drag  Parametrizations 
1  SINGLE- WAVE  AND  SPECTRAL  GRAVITY  WAVE  DRAG  SCHEMES 
An  emerging  trend  in  models  with  middle  atmospheres  is  to 
use  spectral  schemes  to  parametrize  collectively  non-oro- 
graphic  GWD,  with  orographic  GWD  parametrized  separate¬ 
ly  as  a  single-wave  scheme  (e.g.,  Fomichev  et  al.,  2002).  As 
parametrizations  of  non-orographic  sources  such  as  convec¬ 
tion  improve,  it  may  be  natural  to  implement  them  as  single¬ 
wave  schemes  also  (see  Section  3b).  Observationally  at  least, 
this  approach  seems  reasonable:  the  gravity-wave  field  in  the 
middle  atmosphere  shows  “ambient”  wave  spectra  at  most 
locations  (e.g.,  Fritts,  1989;  Eckermann,  1995),  punctuated  by 
intermittent  bursts  of  wave  energy  from  individual  wave 
events  from  sources  such  as  orography  (e.g.,  Eckermann  and 
Preusse,  1999)  and  convection  (e.g.,  Dewan  et  al.,  1998).  The 
larger  issue,  then,  is  how  best  to  combine  these  GWD 
schemes  in  a  model. 

The  simplest  way  is  to  run  them  non-interactively  and  sum 
the  separate  GWD  contributions  from  each.  In  contrast  to  the 
horizontally  localized  column  GWD  from  the  single-wave 
schemes  activated  by  spotty  sources  in  the  lower  atmosphere 
(see  Section  7c),  the  GWD  from  the  spectral  scheme  is  gen¬ 
erally  applied  in  every  model  grid  box  globally  and  thus  has 
a  smoother  overall  structure.  In  such  implementations,  the 
spectral  scheme  is  most  likely  to  be  used  as  a  final  global  tun¬ 
ing  mechanism  for  the  model's  middle  atmosphere,  account¬ 
ing  for  unparametrized  sources  of  GWD  in  much  the  same 
way  that  Rayleigh  friction  is  currently  used  to  stabilize  many 
models  (e.g.,  Jackson  et  al.,  2001):  unlike  Rayleigh  friction, 
however,  well-tuned  spectral  schemes  should  not  lead  to 
unrealistic  reductions  in  model  variability  (e.g.,  Lawrence, 
1997a).  Most  model  tests  to  date  have  implemented  spectral 
and  single-wave  schemes  non-interactively  (e.g.,  McFarlane 
et  al.,  1997;  Manzini  and  McFarlane,  1998;  Scaife  et  al., 
2000;  Osprey,  2001;  Akmaev,  2001;  Scinocca,  2002; 
Fomichev  et  al.,  2002),  with  the  spectral  GWD  scheme  ini¬ 
tialized  and  tuned  using  a  highly  simplified  source  (e.g.,  zon- 
ally  uniform  source-level  spectrum). 

While  this  is  a  sensible  way  to  start,  most  spectral  schemes 
advocate  “interactive  implementations”  in  which  individual 
wave  source  information  is  coupled  within  the  spectral  scheme 
itself  (e.g.,  Hines,  1997b;  Alexander  and  Dunkerton,  1999; 
Medvedev  and  Klaassen,  2000).  Osprey  (2001)  has  tested  both 
interactive  and  non-interactive  implementations  of  orographic 
and  spectral  GWD  schemes  in  offline  single-column  tests  and 
within  a  3D  stratosphere-mesosphere  model.  Differences  in 
each  implementation  were  not  huge,  though  sensitivity  to  the 
strength  of  the  source  was  apparent.  More  studies  along  these 
lines  are  needed  to  clarify  the  dependencies. 

The  distinction  between  the  single-wave  and  spectral 
schemes  reflects,  to  some  extent,  differing  needs  in  various 
models  and  different  theoretical  approaches  to  the  parame- 
trization  problem.  Most  schemes  have  been  implemented 
successfully,  to  some  extent,  in  models  and  there  has,  to  date, 
been  little  effort  to  couple  several  schemes  systematically 
within  a  global  model.  Some  single-wave  GWD  schemes, 


which  were  developed  originally  for  lower  atmospheric  mod¬ 
els,  have  also  been  extended  to  act  as  multiple-wave  (multi¬ 
ple  single-wave)  schemes  to  cover  the  middle  and  upper 
atmospheres  (e.g.,  Jackson,  1993;  Norton  and  Thuburn, 
1999).  Whether  more  explicit  coupling  among  different 
schemes  is  useful  or  even  desirable  remains  to  be  seen:  for 
example,  the  lack  of  success  in  coupling  lower-level  and 
upper-level  radiation  parametrizations  in  models  may  augur 
poorly  for  similar  efforts  with  GWD  schemes. 

2  OROGRAPHIC  GRAVITY  WAVE  DRAG  SCHEME  AND  ENHANCED 
SURFACE  FRICTION  DRAG 

The  surface  friction  drag  induced  by  flows  over/around  com¬ 
plex  orography  has  been  actively  studied:  see,  for  example, 
recent  reviews  by  Belcher  and  Hunt  (1998)  and  Wood  (2000). 
Like  boundary  layers  over  flat  terrain,  studies  suggest  near¬ 
surface  flow  over  hilly  terrain  can  also  be  represented  as  a 
“form  drag”  with  enhanced  effective  roughness  lengths  due  to 
both  vegetative  and  subgrid-scale  orographic  effects  (Wood 
and  Mason,  1993).  In  addition  to  drag  on  the  large-scale  flow, 
surface  friction  can  significantly  affect  orographic  GWD: 
studies  suggest  it  reduces  mountain  wave  amplitudes  and 
breaking  aloft  (Section  7d)  and  is  critical  for  the  realistic  for¬ 
mation  of  downslope  windstorms  and  rotor  flows  (e.g., 
Richard  et  al.,  1989;  Olafsson  and  Bougeault,  1997; 
Leutbecher  and  Volkert,  2000;  Doyle  et  al.,  2000;  Doyle  and 
Durran,  2002).  This  suggests  a  need  to  couple  surface  friction 
and  orographic  GWD  parametrizations  in  models,  an  issue 
that  is  beginning  to  receive  attention  (Wood  et  al.,  2001). 

Inclusion  of  both  a  parametrization  of  an  enhanced  surface 
friction  drag  and  an  orographic  GWD  parametrization  in  large- 
scale  models  can  yield  beneficial  impacts  (Mason,  1986;  Milton 
and  Wilson,  1996;  Gregory  et  al.,  1998;  Hogan  et  al.,  1999).  In 
particular,  Gregory  et  al.  (1998)  showed  that  insertion  of  a  sur¬ 
face  friction  scheme  reduced  parametrized  orographic  gravity 
wave  fluxes  by  half,  possibly  due  to  direct  frictional  retardation 
of  surface  flow  forcing  of  waves.  Such  work  may  help  further 
alleviate  the  excess  of  upper-level  to  lower-level  orographic 
GWD  ratio  noted  in  Section  7d.  These  improvements  follow 
from  running  the  schemes  non-interactively.  Whether  more 
explicit  coupling  is  needed  to  simulate  other  effects  of  surface 
friction  on  orographic  GWD  remains  to  be  seen. 

Enhanced  friction  drag  (i.e.,  form  drag),  however,  should 
be  inserted  into  models  with  caution.  As  discussed  in 
Section  2a,  introduction  of  a  new  drag  mechanism  will  mod¬ 
ify  the  balance  among  the  models’  various  drag  mechanisms. 
For  example,  an  increase  in  surface  friction  drag  may  com¬ 
pensate  for  the  underestimation  of  other  drag  terms,  such  as 
GWD  or  mountain  drag  (see  Kim  and  Hogan,  unpublished 
manuscript,  2003),  leading  to  an  apparent  (i.e.,  illusory) 
improvement  of  the  model  simulation.  Boer  and  Lazare 
(1988)  reported  a  paradoxical  increase  in  net  surface  drag 
(mountain  drag  +  friction  drag)  when  a  GWD  parametrization 
was  removed  from  their  model,  but  this  may  be  due  to  the 
omission  of  GWD  from  their  budget  calculations  (Kim  and 
Hogan,  unpublished  manuscript,  2003). 


Gravity-Wave  Drag  Parametrization  for  Climate  and  NWP  Models  /  83 


f  Interactions  among  Gravity  Wave  Drag  and 
Convection/Cloud  Parametrizations 

As  noted  by  Shepherd  (2000),  subgrid-scale  GWD  schemes 
control  the  accuracy  of  global  middle  atmosphere  simulations 
in  much  the  same  way  that  subgrid-scale  parametrizations  of 
clouds/convection  (and  radiation/land  surface  processes) 
impact  simulations  of  the  global  weather  and  climate  of  the 
troposphere.  It  is  not  surprising,  then,  that  as  global  models 
extend  their  boundaries  to  encompass  both  the  troposphere 
and  middle  atmosphere,  a  number  of  potentially  important 
interactions  among  GWD  schemes  and  convection/cloud 
schemes  are  coming  to  light.  We  consider  a  few  here. 

1  SOURCE  TERMS  FOR  CONVECTIVE  GRAVITY  WAVE  DRAG 
PARAMETRIZATION 

As  discussed  in  Section  3b,  output  from  convective  parame¬ 
trizations  is  beginning  to  be  used  to  initialize  convective 
GWD  parametrizations.  Various  parameters  have  been  used 
to  specify  this  wave  source  term  to  date:  convective  mass 
flux,  precipitation,  latent  heating,  and  so  on.  Since  convec- 
tively  generated  gravity  waves  are  generally  non-stationary, 
formulating  the  radiated  gravity  wave  stress  at  the  cloud  tops 
is  nontrivial.  Steady-state  formulations  of  cloud-top  wave 
stress  induced  by  stationary  diabatic  forcing  by  Chun  and 
Baik  (1998)  are  based  on  a  reference  frame  moving  with  the 
convection,  generating  waves  that  are  stationary  relative  to 
the  diabatic  forcing.  However,  recent  cloud-resolving  numer¬ 
ical  modelling  studies  (e.g.,  Alexander  et  al.,  1995;  Piani  et 
ah,  2000;  Lane  et  ah,  2001)  show  more  momentum  flux  con¬ 
tributed  by  waves  propagating  westward  relative  to  the  con¬ 
vective  system.  These  non-stationary  gravity  waves  need  to 
be  included  in  convective  GWD  parametrizations. 

The  study  of  explicitly  resolved  gravity  waves  generated 
by  tropical  convection  may  aid  this  undertaking.  Such  studies 
in  global  models  (discussed  in  Section  6a)  raised  a  more  gen¬ 
eral  question  about  whether  convective  adjustment  and  cloud 
parametrization  schemes  in  models  are  sufficiently  accurate 
at  present  to  permit  reliable  specifications  of  the  convective 
wave  source  term  in  global  models.  This  is  particularly  criti¬ 
cal  to  assess,  since,  if  unreliable,  the  large  amounts  of 
resolved  wave  activity  generated  by  these  schemes  in  the 
tropics  could  negatively  impact  climate  and  forecast  skill  as 
model  resolutions  increase. 

2  SUBGRID-SCALE  GRAVITY  WAVE  FEEDBACKS  ON  CLOUD  FIELDS 
Fractional  cloud  cover  and  precipitation  within  a  model  grid 
box  are  often  calculated  in  convective  schemes  with  the  aid  of 
probability  distributions  of  water  saturation  due  to  subgrid- 
scale  variability.  While  it  has  generally  been  assumed  that  this 
variability  stems  from  convective  eddies  and  turbulence,  it  is 
clear  that  gravity  waves  and  other  sources  of  subgrid-scale 
variability  at  mesoscale  resolutions  also  contribute  (Balaji 
and  Redelsperger,  1996). 

The  general  banding  of  cloud  cover  by  gravity  waves  is 
well  established  from  observations  (e.g.,  Erickson  and 
Whitney,  1973;  Worthington,  2001),  while  modelling  has 


shown  how  typical  gravity  wave  temperature  variability 
affects  microphysics  of  high  cirrus  (Jensen  and  Toon,  1994) 
and  polar  stratospheric  clouds  (e.g.,  Bacmeister  et  al.,  1999; 
Dornbrack  et  al.,  2001).  Gravity  waves  associated  with  con¬ 
vective  weather  feed  back  to  modify  or  organize  severe  con¬ 
vection  and  precipitation  (e.g.,  Koch  and  Siedlarz,  1999). 

Flows  over  subgrid-scale  orography  also  affect  precipita¬ 
tion.  Smith  (1979)  cites  three  essentially  independent  mecha¬ 
nisms  for  producing  orographic  rain:  large-scale  upslope  rain 
due  to  condensation  caused  by  orographic  ascent  or  oro- 
graphically-triggered  convection;  rainfall  enhancement  by 
low-level  clouds  over  small  hills  through  washout  of 
upper-level  cloud  droplets;  local  downslope  rain  from  cumu¬ 
lonimbus  forming  in  thermals  as  mountain  slopes  are  heated 
by  the  sun.  Some  studies  suggest  orographic  rain  enhance¬ 
ments  are  significantly  impacted  by  mountain  wave  dynamics 
(Richard  et  al.,  1987;  Reinking  et  al.,  2000;  Brady  and 
Waldstreicher,  2001). 

Recent  work  has  focused  on  the  way  in  which  gravity 
waves  and  other  sources  of  unresolved  mesoscale  variability 
can  significantly  influence  output  of  cloud/convection  param¬ 
etrizations  in  large-scale  models  (Balaji  and  Redelsperger, 
1996;  Cusack  et  al.,  1999;  Nilsson  et  al.,  2000).  In  models  it 
is  conceivable  that  GWD  schemes  could  be  used  to  feed  sub- 
grid-scale  gravity  wave  fluctuations  directly  into  cloud 
parametrizations  (e.g.,  Bacmeister  et  al.,  1999). 

3  DIABATIC  CLOUD/MOISTURE  EFFECTS  ON  OROGRAPHIC 
GRAVITY  WAVE  DRAG 

Diabatic  effects  within  clouds,  such  as  latent  heat  release  or 
cloud-top  infra-red  cooling,  can  modify  the  tropospheric 
environment  in  ways  that  significantly  impact  the  trapping 
and  amplitude  evolution  of  mountain  waves  and  the  GWD 
they  produce  (e.g.,  Durran  and  Klemp,  1982a,  1982b; 
Weissbluth  and  Cotton,  1989;  Doyle  and  Shapiro,  2000; 
Doyle  and  Smith,  2003).  Surgi  (1989)  included  a  moisture- 
modified  Brunt- Vaisala  frequency  tuned  arbitrarily  to  various 
cloud  types  in  an  orographic  GWD  scheme  that  was  tested  in 
a  global  model.  Since  then,  these  effects  have  received  little 
attention  in  parametrizations.  Again,  the  issue  of  whether 
convective  schemes  in  models  are  good  enough  at  present 
(Section  6a)  may  dictate  the  feasibility,  or  otherwise,  of  pur¬ 
suing  developments  along  these  lines. 

g  Interactions  between  Resolved  Gravity  Waves  and 
Parametrized  Gravity  Wave  Drag  in  High-Resolution 
Models 

As  discussed  in  Section  6,  increases  in  resolution  are  allow¬ 
ing  models  to  resolve  the  “outer  scales”  of  the  gravity  wave 
spectrum  and  to  capture  some  GWD.  This  implies  that 
any  changes  in  grid-box  resolution  in  next-generation  models 
will  be  accompanied  by  changes  in  resolved  GWD  that 
will  necessitate  retuning  of  the  parametrized  subgrid-scale 
GWD.  Beyond  such  practical  concerns,  this  raises  some 
deeper  and  perhaps  thornier  issues  for  subgrid-scale  GWD 
parametrization. 
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1  TEMPORAL  VARIABILITY  OF  GRAVITY  WAVES 

At  each  model  time  step,  column  parametrizations  of  GWD 
ingest  a  vertical  profile  of  model  winds  and  temperatures 
within  a  given  grid  box  and  propagate  gravity  waves  through 
this  vertical  column  as  if  the  atmosphere  were  static.  Of 
course  these  winds  and  temperatures  vary  over  time,  but  if 
wave  propagation  times  are  sufficiently  fast  this  complication 
can  usually  be  ignored. 

As  models’  spatial  resolutions  improve,  however,  model 
time  steps  Ar  decrease  unless  a  more  efficient  transport  algo¬ 
rithm  such  as  a  semi-Lagrangian  scheme  is  used,  and  thus 
there  is  greater  potential  for  significant  temporal  atmospheric 
variability  to  be  resolved.  Furthermore,  as  the  tops  of  models 
are  extended  into  the  middle  atmosphere  and  beyond,  the  total 
time  required  for  any  given  gravity  wave  to  propagate  verti¬ 
cally  through  the  column  from  its  tropospheric  source,  tpROp, 
increases  considerably.  Decreases  in  At  and  increases  in  tpROp 
both  conspire  to  make  the  neglect  of  temporal  variability  in 
modelling  gravity  wave  propagation  in  GWD  parametriza¬ 
tions  potentially  less  tenable.  In  situations  where  time  depen¬ 
dence  of  the  background  atmosphere  has  been  shown  to  be 
important,  the  modifications  it  makes  to  GWD  profiles  can  be 
quite  large  (see,  e.g.,  Eckermann  and  Marks,  1996). 

Increased  temporal  variability  in  models  also  leads  to 
greater  time  variations  in  applied  subgrid-scale  GWD.  Steady 
GWD  drives  the  atmosphere  to  a  new  steady  dynamical  state, 
whereas  temporally  unsteady  GWD  produces  vacillating 
mean  flows.  In  particular,  if  steady  GWD  in  some  region  sud¬ 
denly  increases  or  abates,  the  flow  there  becomes  unbalanced 
and  will  try  to  relax  to  a  new  balanced  state  through  reradia¬ 
tion  of  secondary  gravity  waves  (e.g.,  Vadas  and  Fritts,  2001). 
It  could  be  argued  that  spectral  models  parametrize  the  net 
effects  of  all  waves,  including  reradiated  secondary  waves, 
and  as  such  might  also  be  less  affected  by  the  “transit  time” 
issue  for  single  waves  discussed  here.  As  model  resolutions 
improve,  temporal  intermittency  in  GWD  applied  by  subgrid- 
scale  parametrizations  could  cause  these  models  to  radiate  an 
explicitly  resolved  spectrum  of  secondary  gravity  waves 
spontaneously  much  of  the  time.  Their  existence  and  any 
potential  effects  on  model  accuracy  are  unclear  at  present. 

2  WAVE-WAVE  INTERACTIONS 

Much  of  the  increased  spatio-temporal  variability  in  high-res- 
olution  models  is  produced  by  explicitly  resolved  long  wave¬ 
length  gravity  waves,  particularly  above  the  tropopause 
(Koshyk  and  Hamilton,  2001;  Section  7).  These  resolved 
waves  produce  oscillations  or  “wiggles”  on  the  background 
wind  and  temperature  profiles  passed  to  GWD  parametriza¬ 
tions.  In  models  where  resolved  waves  are  energetic, 
these  oscillations  can  significantly  modulate  the  predicted 
subgrid-scale  GWD,  and  thus  the  GWD  scheme  implicitly 
simulates  some  form  of  wave-wave  interaction  between  the 
resolved  gravity  waves  and  the  subgrid-scale  (parametrized) 
gravity  waves. 

Such  situations  were  never  envisaged  in  the  formulation  of 
most  GWD  parametrizations,  and  thus  these  effects  require 


scrutiny.  Detailed  modelling  of  the  modulations  of  short  grav¬ 
ity  waves  propagating  through  longer  gravity  waves  indicate 
that  the  interactions  are  much  more  complicated  than  single¬ 
wave  GWD  parametrizations  would  predict:  for  instance,  the 
time  dependence  of  the  long  wave  field  usually  cannot  be 
neglected  when  modelling  the  propagation  and  amplitudes  of 
short  waves  (Eckermann  and  Marks,  1996;  Walterscheid, 
2000).  Indeed,  a  whole  class  of  potentially  relevant  wave- 
wave  interactions  may  occur,  none  of  which  is  included  in 
single-wave  GWD  schemes.  While  some  spectral  models  are 
formulated  to  encapsulate  the  net  spectral  effects  of  certain 
wave-wave  interactions,  the  GWD  parametrizations  they 
spawn  are  again  formulated  under  the  premise  that  they  act 
on,  and  respond  to,  only  a  large-scale  atmosphere  devoid  of 
resolved  gravity  waves  (e.g.,  Hines,  1997a,  1997b;  Medvedev 
and  Klaassen,  2000). 

In  short,  existing  GWD  parametrizations  are  not  designed 
for  models  containing  energetic  resolved  gravity  waves.  As 
more  pressing  parametrization  problems  are  tackled  and 
solved,  this  wave-wave  interaction  issue  may  become  pro¬ 
gressively  more  important,  particularly  as  model  resolutions 
improve  and  upper  boundaries  are  extended.  One  near-term 
solution  could  follow  from  modifying  Lander  and  Hoskins’ 
(1997)  suggestion  to  use  only  the  “believable  scales”  of  a 
model  as  input  for  parametrizations.  Current  GWD  parame¬ 
trizations  in  high-resolution  models  could  remain  “believ¬ 
able”  if  they  used  vertical  profiles  of  model  winds  and 
temperatures  that  have  been  “averaged"  or  “filtered"  so  as  to 
eliminate  resolved  gravity  wave  oscillations.  Studies  of 
resolved  wave-wave  interactions  within  models  could  also  be 
very  helpful  in  understanding  and  ultimately  parametrizing 
these  effects. 

h  The  Jet  Stream  ( as  an  Example  of  Other  Gravity  Wave 
Sources ) 

Gravity  wave  sources  besides  convection  and  topography 
exist,  which  have  proven  difficult  to  understand  and  parame¬ 
trize.  Section  3  of  Fritts  and  Alexander  (2003)  provides  a 
good  review  of  a  variety  of  sources  that  merit  attention.  In  this 
subsection,  we  discuss  waves  from  the  “jet  stream”,  since 
their  potential  relevance  to  GWD  in  the  summer  middle 
atmosphere,  for  instance,  has  been  established  (e.g.,  Btihler 
and  McIntyre,  1999;  Scinocca  and  Ford,  2000;  see  also 
Fig.  1).  In  the  lower  atmosphere,  they  may  also  exert  an 
important  “emission  drag”  on  the  jet  itself  (Sutherland  and 
Peltier,  1995). 

While  observational  associations  go  back  many  years  (e.g., 
Hooke  and  Hardy,  1975),  modelling  jet-stream  gravity  waves 
has  proven  difficult.  Uccellini  and  Koch  (1987)  tackled  the 
problem  from  an  observational  meteorological  perspective. 
They  noticed  that  energetic  “mesoscale  waves”  reported  over 
the  continental  U.S.A.  arose  in  conjunction  with  diffluent 
“exit”  regions  of  the  jet  at  -300  hPa.  They  argued  that 
geostrophic  adjustment  of  these  regions  radiates  long-wave- 
length  inertia  gravity  waves.  In  their  theory,  these  waves 
become  trapped  in  the  troposphere  due  to  a  critical  layer  made 
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reflective  (or  “overreflective”)  by  unstable  vertical  shear. 
Subsequent  observations  and  modelling  have  supported  the 
general  scenario,  while  substantially  developing  the  overall 
picture,  and  with  details  varying  on  a  case-by-case  basis  (e.g., 
Koch  and  Dorian,  1988;  Powers,  1997;  Koch  and  Siedlarz, 
1999;  Zhang  et  al.,  2001).  Their  potential  GWD  contribution 
is  unclear.  In  the  absence  of  ducting  by  an  unstable  critical 
layer,  these  waves  should  propagate  into  the  middle  atmos¬ 
phere,  as  has  been  verified  by  observation  (Guest  et  al.,  2000; 
Pavelin  et  al.,  2001)  and  3D  modelling  studies  (Luo  and  Fritts, 
1993;  O’Sullivan  and  Dunkerton,  1995).  Breaking  of  these 
waves  higher  up  (e.g.,  Pavelin  et  al.,  2001)  produces  GWD. 

However,  other  types  of  instabilities  can  also  radiate  waves 
from  the  jet;  for  example,  the  unstable  shear  layers  that  duct 
waves  in  the  Uccelini  and  Koch  (1987)  model  can  also  radi¬ 
ate  waves  via  dynamical  (Kelvin-Helmholtz)  instabilities  in 
certain  circumstances  (e.g.,  Sutherland  and  Peltier,  1995; 
Btihler  et  al.,  1999;  Scinocca  and  Ford,  2000).  Still  other 
models  associate  jet-level  wave  generation  with  other  dynam¬ 
ics,  such  as  “inertial  instability”  (Ciesielski  et  al.,  1989)  or 
“Lighthill  radiation”  (Ford,  1994;  Medvedev  and  Gavrilov, 
1995;  Reeder  and  Griffiths,  1996;  Ford  et  al.,  2000). 

These  generation  mechanisms  are  all  interrelated  to  varying 
extents,  focusing  on  some  diagnosed  form  of  dynamical 
imbalance  or  instability  that  develops  in  the  jet  and  then  trig¬ 
gers  radiation  of  gravity  waves  to  restore  balance.  A  unified 
theory  of  such  dynamics  has  been  sought  but  remains  elusive; 
for  example,  the  formal  existence  or  otherwise  of  a  “slow  (or 
quasi-slow)  manifold”  that  would  unambiguously  separate 
“balanced”  motion  from  “unbalanced”  motions  such  as  grav¬ 
ity  waves  (e.g.,  McIntyre  and  Norton,  2000;  Ford  et  al.,  2000). 
A  variety  of  measures  have  been  proposed  for  diagnosing  bal¬ 
ance  to  various  levels  of  approximation  for  potential  use  in 
identifying  when  and  where  waves  might  originate  due  to  jet 
imbalances  (e.g.,  Knox,  1997;  Olsson  and  Cotton,  1997; 
Zhang  et  al.,  2000).  Koch  and  O’Handley  (1997)  showed  that 
some  of  these  measures,  when  calculated  from  NWP  model 
output,  showed  promise  in  forecasting  mesoscale  wave  occur¬ 
rences  in  the  troposphere.  Other  quantities  (e.g.,  Rossby  radii 
of  deformation,  duct  depth)  seem  to  specify  observed  para¬ 
meters  of  these  waves  reasonably  well  (Uccellini  and  Koch, 
1987;  Koch  and  O’Handley,  1997).  Such  terms  could  be  cal¬ 
culated  within  global  models  as  source  terms  for  a  subgrid- 
scale  parametrization  of  “jet-exit”  GWD. 

One  difficulty,  however,  may  be  that,  similar  to  wave- 
induced  instabilities  (see  Section  6c2),  these  source  regions  of 
flow  instability  may  often  be  of  subgrid  scales  themselves 
(see  e.g.,  Sutherland  and  Peltier,  1995;  Section  5b  of  Scinocca 
and  Ford,  2000),  making  this  source  of  GWD  hard  to  param¬ 
etrize.  If  so,  the  only  recourse  may  be  higher  model  resolu¬ 
tions.  O’Sullivan  and  Dunkerton  (1995)  used  a  simplified 
global  spectral  model  to  simulate  gravity  wave  radiation 
explicitly  from  exit  regions  of  evolving  mid-latitude  baroclin- 
ic  disturbances.  Appreciable  wave  radiation  occurred  at  T106 
with  Az-700  m  in  the  stratosphere.  Comparable  or  better  res¬ 
olution  is  used  in  many  NWP  models,  which  suggests  (see 


Section  6)  that  these  jet-stream  source  regions  and  some  of 
the  large-scale  gravity  waves  they  radiate  may  already  be 
appearing  explicitly  in  some  higher  resolution  global  models. 

8  Conclusions 

Expanding  scientific  knowledge  and  advancing  computing 
technologies  have  encouraged  converging  trends  in  develop¬ 
ment  of  global  atmospheric  models.  For  example,  the  distinc¬ 
tions  between  global  NWP  models  and  GCMs  are  now 
becoming  less  obvious:  NWP  models  are  gradually  increasing 
the  temporal  range  of  their  numerical  integrations,  while 
GCMs  are  gradually  increasing  their  spatial  resolutions. 
Further,  middle  atmosphere  modellers  have  come  to  appreci¬ 
ate  the  importance  of  a  realistic  troposphere  (e.g.,  Hamilton, 
1996;  Garcia,  2000),  while  tropospheric  modellers  and  fore¬ 
casters  are  incorporating  middle  atmospheres  to  improve  both 
climate  and  forecast  skill  (e.g.,  Pawson  et  al.,  2000;  Thompson 
et  al.,  2002).  As  spatial  model  resolutions  increase  the  need  for 
non-hydrostatic  model  formulations  has  arisen  (see  Fig.  11; 
see  also  Smolarkiewicz  et  al.,  2001),  which  removes  a  major 
distinction  between  global  and  regional  NWP  models  (e.g., 
UK  Met  Office  Unified  Model;  Cullen  et  al.,  1997). 

There  will  still  be  important  differences  of  course;  for 
example,  the  tops  of  NWP  models  are  expected  to  be  extend¬ 
ed  to  include  the  middle  atmosphere,  but  perhaps  not  as  high 
as  those  of  GCMs  which  include  the  upper  atmosphere.  NWP 
models  will  always  tend  to  have  higher  resolutions  and  short¬ 
er  integration  time  steps  than  GCMs,  with  GCMs  including 
more  diverse  physics  and  chemistry  packages  and  NWP  mod¬ 
els  focusing  more  on  data  assimilation. 

Nonetheless,  these  evolving  similarities  mean  that  differing 
perspectives  on  GWD  parametrization  that  have  emerged 
among  various  segments  of  the  global  modelling  communi¬ 
ties  (e.g.,  troposphere  versus  middle  atmosphere  modellers; 
Sections  3  and  4)  now  need  to  be  combined  and  reconciled  for 
wider  future  modelling  use.  This  paper  has  attempted  to  sum¬ 
marize  some  of  these  perspectives  and  some  of  the  issues  that 
arise  in  combining  them  more  generally.  Emerging  issues 
have  also  been  discussed  and/or  anticipated,  such  as  explicit¬ 
ly  resolved  gravity  waves  and  coupling  among  parametriza- 
tions  (Sections  6  and  7).  We  also  speculated  on  the  model 
resolutions  ultimately  required  to  eliminate  the  need  for 
GWD  parametrization  (Section  6). 

The  science  of  GWD  parametrization  is  only  about  20 
years  old,  yet  it  has  progressed  and  evolved  significantly  dur¬ 
ing  that  time.  Thus,  beyond  the  general  issues  and  trends  cited 
above,  it  is  difficult,  and  probably  futile,  to  speculate  on 
exactly  how  things  might  develop  in  the  next  5-10  years  as 
next-generation  global  models  come  online.  It  is  certainly 
safer  to  assert  that  many  of  the  current  uncertainties  in  GWD 
parametrization  specifically,  and  atmospheric  dynamics  gen¬ 
erally,  may  be  solved  in  no  small  measure  by  the  insights  pro¬ 
vided  by  these  next-generation  modelling  tools.  Such 
breakthroughs  imply  much  improved  understanding  and  fore¬ 
casting  of  our  atmosphere  as  a  whole,  and  in  particular  the 
role  of  gravity  wave  processes  within  it. 
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Appendix  A:  Basic  theories  and  nomenclature  relevant 
to  gravity  wave  drag  parametrization 

We  will  briefly  summarize  some  basic  concepts  and  dynami¬ 
cal  theories  that  are  needed  to  understand  the  parametrization 
of  GWD  in  large-scale  atmospheric  models.  Figure  10  pro¬ 
vides  an  idealized  picture  of  the  type  of  gravity  wave  genera¬ 
tion,  propagation  and  dissipation  processes  in  the  atmosphere 
that  parametrizations  seek  to  approximate.  Parametrizations 
focus  on  the  drag  effects  produced  by  gravity  wave  dissipa¬ 
tion  that  are  important  for  NWP  and  climate  models.  This  rep¬ 
resents  the  end  stage  of  the  wave’s  life  cycle,  but  to  describe 
it  accurately,  parametrizations  must  also  describe  the  earlier 
“birth”  and  “evolution”  stages  that  lead  up  to  it.  This  basic 
three-stage  “life  cycle”  paradigm  for  a  gravity  wave  under¬ 
pins  the  construction  of  most  GWD  parametrization  schemes 
(see,  e.g..  Section  7  of  Fritts  and  Alexander,  2003),  and  will 
be  used  here  to  introduce  common  terminology  and  concepts 
encountered  in  the  field. 


a  Surface  Drag 

We  begin  by  briefly  discussing  drag  processes  operating  at 
the  surface.  Even  for  stably  stratified  flow  over  a  symmetric 
frictionless  surface  obstacle,  a  pressure  difference  can  occur 
between  the  surface  flow  upstream  and  downstream  of 
an  obstacle.  For  a  2D  ridge,  this  yields  a  “drag”  force 


x  = 


r  dh{x) 

)  Psurf\x)  , 


dx  where  Pulrf(x)  is  the  surface  pressure 


as  a  function  of  horizontal  distance  x  in  response  to  variations 
in  elevation  due  to  the  ridge  height  h(x),  and  the  integral  is 
across  the  full  x  extent  of  the  obstacle.  This  drag  force, 
referred  to  as  “mountain  drag”  or  sometimes  “pressure  drag” 
(e.g..  Smith,  1978),  represents  a  transfer  of  angular  momen¬ 
tum  from  the  atmosphere  to  the  solid  Earth,  which  acts  to 
decelerate  the  zonal  flow  at  the  surface. 


For  higher  mountain  elevations  h(x)  the  flow  sometimes 
cannot  pass  over  the  mountain,  and  thus  the  flow  becomes 
blocked  upstream  for  2D  ridges  or  flows  around  the  mountain 
peak  for  3D  obstacles,  leading  to  drag  and  altered  atmospher¬ 
ic  responses  (e.g.,  lee  vortices). 

There  must  be  an  equal  and  opposite  exchange  of  momen¬ 
tum  from  the  solid  Earth  back  to  the  atmosphere.  This  atmos¬ 
pheric  response  manifests  itself  in  different  forms  depending 
on  the  wind  speed  and  horizontal  and  vertical  scales  of  the 
orographic  feature  (e.g.,  Smith,  1978).  At  planetary  topo¬ 
graphic  scales,  for  example,  the  response  can  be  in  the  form 
of  radiated  planetary  Rossby  waves,  a  process  that  is  resolved 
in  most  GCMs  and  NWP  models.  At  very  short  topographic 
scales  the  atmospheric  response  is  turbulent,  often  parame¬ 
trized  in  terms  of  an  effective  surface  roughness  that  acts  like 
a  frictional  force  on  the  flow  within  the  planetary  boundary 
layer  (“form  drag”).  These  processes  are  discussed  in  more 
detail  in  Sections  7b,  7d  and  7e2  (see  also  recent  reviews  by 
Belcher  and  Hunt  (1998),  Wood  (2000)  and  Smith  (2001)).  At 
horizontal  topographic  scales,  -10-100  km,  where  buoyancy 
forces  tend  to  dominate,  the  atmospheric  response  is  often  a 
radiated  field  of  gravity  waves  (mountain  waves)  with  a  ver¬ 
tical  flux  of  horizontal  momentum  density  pu'w'  that  bal¬ 
ances  the  surface  mountain  drag  x  produced  by  these 
particular  topographic  features,  where  p  is  the  atmospheric 
density,  u  and  w'  are,  respectively,  the  horizontal  and  vertical 
velocity  perturbations  of  the  wave,  and  the  overbar  denotes 
averaging  over  a  full  wave  cycle  (wave  momentum  flux  is 
discussed  further  in  Section  dl  of  this  Appendix).  Since 
10-100-km  scale  topography  is  not  fully  resolved  in  GCMs  or 
NWP  models,  this  subgrid-scale  surface  mountain  drag  and 
the  radiated  mountain  waves  must  be  parametrized. 

b  Wave  Generation 

1  OROGRAPHIC  VERSUS  NON-OROGRAPHIC  GRAVITY  WAVES 

As  we  have  just  noted,  one  very  important  source  of  gravity 
waves  for  drag  parametrization  at  all  levels  is  flow  over 
topography  (“orographic”  forcing),  and  drag  due  to  mountain 
waves  has  been  widely  parametrized  (details  discussed  in 
Section  3a).  In  fact,  gravity  waves  are  generated  by  a  variety 
of  other  sources  that  are  located  mostly  in  the  troposphere, 
however,  waves  from  sources  other  than  mountains  are  less 
well  understood,  and  currently  are  often  ignored  or  parame¬ 
trized  collectively.  For  this  reason,  GWD  parametrizations 
are  often  classified  as  either  “orographic”  or  “non-orograph- 
ic”.  Sources  of  non-orographic  gravity  waves  include  con¬ 
vection  (Sections  3b,  6a),  jet  streams  (Section  7h)  and  fronts, 
among  others.  See  Section  3.1  of  Fritts  and  Alexander  (2003) 
for  a  detailed  review  of  atmospheric  processes  that  generate 
gravity  waves. 

2  SINGLE  GRAVITY  WAVES  VERSUS  WAVE  SPECTRA 
Parametrizations  tend  to  fall  into  “single-wave”  schemes, 
which  deal  with  individual  waves,  and  “spectral”  schemes, 
which  collectively  parametrize  the  effects  of  a  broad  spec¬ 
trum  of  gravity  waves.  Based  on  an  earlier  understanding  that 
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Fig.  10  Schematic  of  various  gravity  wave  production,  propagation  and  dissipation  processes  that  parametrizations  seek  to  capture,  along  with  associated 
processes  (see  symbol  key  at  top).  Surface  orography,  convection  and  jet  streams  are  shown  as  sources  of  waves,  with  stationary  mountain  waves 
launched  quasi-vertically  and  non- stationary  waves  launched  obliquely  into  the  atmosphere.  Wave  breaking  is  depicted  along  wave  propagation  paths 
in  various  regions  of  the  atmosphere.  Wave  breaking  is  more  prevalent  higher  in  the  atmosphere,  where  wave  propagation  paths  cross  and  thus  waves 
interact  more,  yielding  broader  wave  spectra,  and  some  breaking  reradiates  secondary  waves.  Other  waves  break  lower  down,  nearer  their  source  or 
reflect  vertically  and  break  via  trapping/resonance  effects  (see  mountain  wave,  lower-right).  The  figure  depicts  a  basic  paradigm  of  the  troposphere  as 
a  source  region,  the  stratosphere  as  a  propagation  region  and  the  upper  stratosphere  and  mesosphere  as  a  strong  breaking  region. 
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gravity  wave  sources  mostly  radiate  coherent  wave  packets 
(Fig.  10),  most  near-source  parametrizations  have  been  sin¬ 
gle-wave  schemes  (though  not  always  -  see,  e.g.,  Gregory  et 
al.,  1998).  As  gravity  waves  propagate  away  from  their  tro¬ 
pospheric  source  regions,  their  propagation  paths  start  to 
meander  further  and  intersect  with  those  of  other  waves 
(Fig.  10),  so  that  a  broader  distribution  of  wavelengths  and 
frequencies  is  encountered  higher  up  and  total  wave  variance 
is  less  identifiable  with  low-level  sources.  For  these  historical 
and  scientific  reasons,  spectral  GWD  schemes  tend  to  be 
applied  in  the  middle  atmosphere,  whereas  in  the  lower 
atmosphere  single-wave  schemes  tend  to  dominate  GWD 
parametrizations.  Flowever,  these  distinctions  are  oversimpli¬ 
fied  and  somewhat  arbitrary:  for  example,  broad  gravity  wave 
spectra  are  often  observed  in  the  troposphere,  while  coherent 
single-wave  patterns  from  identifiable  sources  sometimes 
extend  well  into  the  middle  atmosphere  (e.g.,  Dewan  et  al., 
1998).  Also,  the  single-wave  schemes  are  not  restricted  to 
treating  just  one  wave,  instead,  it  is  common  to  apply  single¬ 
wave  schemes  iteratively  to  a  set  of  waves,  summing  the  sin¬ 
gle-wave  contributions  from  each  to  yield  the  total  drag  (e.g., 
Norton  and  Thuburn,  1999).  These  and  other  issues  are  dis¬ 
cussed  in  greater  detail  in  Sections  3  and  4. 

3  STATIONARY  VERSUS  NON- STATIONARY  GRAVITY  WAVES 
Source  parametrizations  are  important  because  they  set  the 
properties  of  the  waves  (e.g.,  wavelength,  amplitude)  that 
control  subsequent  wave  evolution  and  the  form  and  intensi¬ 
ty  of  the  wave-induced  drag.  A  particularly  important  para¬ 
meter  for  GWD  calculations  is  the  wave’s  ground-based 
horizontal  phase  speed,  c,  which  is  set  at  the  source  and  typi¬ 
cally  assumed  to  remain  constant  thereafter.  As  discussed  in 
the  next  two  subsections,  its  value  determines  where  waves 
tend  to  break  and  the  strength  of  the  GWD.  Parametrizations 
commonly  differentiate  between  “stationary”  (c  =  0)  and 
“non- stationary”  (c  =£  0)  gravity  waves,  mainly  because  this 
classification  in  parametrizations  is  currently  synonymous 
with  orographic  and  non-orographic  waves,  respectively. 
While  this  correspondence  is  a  useful  guide  for  parametriza- 
tion,  in  reality  sources  other  than  mountains  can  generate  sta¬ 
tionary  waves  (discussed  in  Sections  3b  and  7h),  while 
orography  can  generate  both  stationary  and  non-stationary 
gravity  waves  (e.g.,  Nance  and  Durran,  1997,  1998). 

C  Wave  Propagation  and  Evolution 
1  INTERNAL  VERSUS  EXTERNAL  GRAVITY  WAVES 
The  linearized  2D  equation  governing  vertical  velocities  due 
to  gravity  waves  (neglecting  the  Coriolis  force,  curvature  of 
the  horizontal  wind  U(z)  and  density  effects)  can  be  written  as 

^  +  (e2  -  k2)w  =  o.  (i) 

Here,  w(k,z,t )  is  the  Fourier-transformed  vertical  velocity 
w(x,z,t)  at  height  z,  and  €2(z)  =  N2(z)/[c  -  t/(z)]2,  where 


N(z)  is  the  Brunt-Vaisala  frequency,  and  k  =  2nlXx  is  the 
horizontal  wavenumber.  For  stationary  mountain  waves 
( c  =  0),  f2(z)  is  known  as  the  Scorer  parameter  (e.g.,  Durran, 
1990).  For  an  assumed  form  of  solution  w(k,z,t )  =  W(k) 
exp [i(\m(z)dz  -  kct)],  Eq.  (1)  yields  the  local  gravity  wave  dis¬ 
persion  relation 

m2(z)  =  i\z)  -k2  =  k2[N2/&r(z)  -  1],  (2) 

which  relates  k  to  the  vertical  wavenumber  m  and  intrinsic 
frequency  cb  =  k(c  -  U).  The  form  of  the  gravity  wave 
depends  on  the  sign  of  m2(z)  =  €2(z)  -  k2,  which  is  determined 
mostly  by  flow  conditions  and  the  source  of  the  gravity  wave. 
The  wave  is  propagating  or  “internal”,  if  m 2  =  €2-k2>  0,  i.e., 
m  is  real.  When  nr  =  €2  -k2  <  0  the  solution  is  “external”,  in 
that  m  becomes  imaginary,  the  wave  is  not  propagating,  and 
w(k,z,t)  decays  with  height  as  exp[-flm(z)lc/z].  When  a  free 
propagating  wave  ( m 2  >  0)  propagates  to  a  height  z  where 
m2(zt)— >0,  its  vertical  wavelength  locally  becomes  infinitely 
long  and  wave  energy  is  reflected  downwards,  with  wave 
solutions  becoming  external  above  this  “turning  point”  zf.  In 
the  current  approximation  Eq.  (2),  turning  points  occur  where 
i2  =  k2  or  equivalently  where  cb  =  N2. 

For  stationary  mountain  waves,  turning  points  can  often 
lead  to  the  trapping  of  wave  energy  between  the  surface  and 
the  turning  point.  Figure  1 1  shows  examples  of  these  regimes, 
plotted  as  a  function  of  zonal  wavelength  X  and  height  z 
under  typical  mid-latitude  winter  conditions  for  U(z)  and 
N(z).  The  internal  mountain  waves  to  the  right  of  this  figure 
propagate  obliquely  away  from  their  surface  topography.  The 
waves  with  longer  horizontal  wavelengths  X  propagate  into 
the  stratosphere,  whereas  the  shorter  scale  waves  encounter 
turning  points  zt,  denoted  by  the  thick  €2  =  k2  line.  These 
waves  reflect  vertically  and  become  “trapped”  between  the 
surface  and  z  and  can  be  channelled  downstream  within  their 
duct  layer.  The  magnitude  of  this  trapping  can  be  measured 
by  the  vertical  decrease  of  the  Scorer  parameter,  €2(z). 

2  HYDROSTATIC  VERSUS  NON-HYDROSTATIC  GRAVITY  WAVES 
Gravity  waves  with  £2/k2  »  1  are  termed  “hydrostatic”,  since 
simplified  hydrostatic  equations  adequately  describe  their 
propagation  and  evolution  (e.g.,  Lindzen,  1981).  For  exam¬ 
ple,  the  dispersion  relation  Eq.  (2)  simplifies  to  nr  ~  €2  = 
N2/(c  -  U)2  for  hydrostatic  waves.  Many  parametrization 
schemes  are  based  on  a  hydrostatic  framework.  Conversely, 
internal  gravity  waves  for  which  £2/k2  ~  O(l)  are  termed 
“non-hydrostatic”.  From  this  limit  and  the  dispersion  relation 
Eq.  (2),  we  note  that  reflected  and  ducted  waves  are  non¬ 
hydrostatic,  and  that  free -propagating  non-hydrostatic  waves 
tend  to  have  relatively  shorter  horizontal  wavelengths  X  , 
longer  vertical  wavelengths  X,,  faster  intrinsic  phase  speeds 
Ic  (=  c  -  U) I,  higher  intrinsic  frequencies  cb  compared  to 
hydrostatic  waves  (see  Fig.  11).  Non-hydrostatic  waves  are 
important  in  the  parametrization  of  low-level  “trapped” 
mountain  waves  and  “fast”  gravity  waves  generated  by  con¬ 
vection  (see  Section  3). 
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Fig.  11  Spectrum  of  the  vertical  propagation  of  orographic  gravity  waves  as  a  function  of  height  and  horizontal  scale,  under  typical  mid-latitude  northern  hemi¬ 
sphere  winter  conditions.  €2  denotes  the  Scorer  parameter,  which  depends  on  static  stability,  wind  speed  and  curvature,  and  k  =  2 tUX  denotes  the  zonal 
wavenumber.  (Based  on  Iwasaki  et  al.  (1989;  with  permission  from  the  Meteorological  Society  of  Japan)  and  Kim  (1996;  with  permission  from  the 
American  Meteorological  Society)) 


Inclusion  of  the  Coriolis  parameter  (inertial  frequency) /in 
Eq.  (1)  leads  to  a  dispersion  relation  in  which  the  N2/  b)2  -  I 
factor  in  Eq.  (2)  becomes  ( N 2  - &>2)/(&> 2  -/2).  Thus, /imposes 
a  low-frequency  limit  for  gravity  waves,  so  that  allowable  fre¬ 
quencies  for  internal  gravity  waves  are  1/1  <  Icbl  <  N.  Gravity 
waves  with  frequencies  near  1/1  are  termed  "inertia  gravity 
waves”,  since  the  Coriolis  force  modifies  their  propagation 
and  polarization  characteristics  compared  to  gravity  waves 
for  which  Icbl  »  I/I. 

3  THE  CRITICAL  LEVEL  Zc 

An  important  physical  process  that  needs  to  be  considered  in 
GWD  parametrizations  is  the  “critical  level”,  an  altitude 
where  the  wind  speed  component  along  the  wave’s  horizontal 
wave  vector  equals  its  phase  speed  c:  in  2D  and  ignoring  rota¬ 
tion  the  criterion  is  c  =  U(zc),  where  zc  is  the  critical  level 
(Lindzen,  1981).  Since  c  =  c  -  U(z)  in  this  case,  this  corre¬ 
sponds  to  vanishing  intrinsic  phase  speed  Icl  and  vertical 
wavelength  X_  from  Eq.  (2).  With  rotation  added,  the  equiva¬ 


lent  criterion  is  Icbl  -  |/l  and  Icl  becomes  very  small  rather  than 
vanishing  completely.  In  most  circumstances  gravity  waves 
dissipate  strongly  as  they  approach  critical  levels  and  are 
totally  dissipated  at  or  just  below  the  critical  level.  This 
process  is  important  as  a  wave  filtering  mechanism  for  the 
middle  atmosphere  (see  Section  4a).  However,  strong  wave 
dissipation  and/or  mean  wind  shear  near  a  critical  level  can 
induce  strong  non-linearity  that  can  sometimes  turn  these 
regions  into  efficient  resonant  reflectors  that  vertically  trap 
wave  energy  (e.g.,  Peltier  and  Clark,  1979).  Such  effects  are 
important  in  parametrizing  GWD  due  to  low-level  mountain 
waves  (Section  3a)  and  tropospheric  gravity  waves  associat¬ 
ed  with  the  jet  stream  (Section  7h).  The  effects  of  non-linear 
wave  trapping  can  be  represented  by  the  vertical  derivative  of 
the  Scorer  parameter,  similar  to  the  linear  trapping  case 
(Section  cl  of  this  Appendix).  The  Scorer  parameter  has  been 
used  explicitly  in  some  parametrization  schemes  to  represent 
the  effects  of  linear  and/or  non-linear  wave  trapping  (Kim  and 
Mahrt,  1992;  Kim  and  Arakawa,  1995;  Gregory  et  al.,  1998). 
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d  Wave  Dissipation  and  Gravity  Wave  Drag 
1  MOMENTUM  FLUX  /  E-P  FLUX  DIVERGENCE 
Steady  2D  gravity  waves  conserve  their  vertical  flux  of  hori¬ 
zontal  wave  momentum  flux  density  in  the  absence  of  dissi¬ 
pation.  For  gravity  waves  propagating  east-west,  this  is  given 
by  p  u'w '  where  u  and  vv'  are,  respectively,  the  zonal  and  ver¬ 
tical  velocity  perturbations  of  the  wave.  Eliassen  and  Palm 
(1961)  first  derived  this  for  orographic  gravity  waves,  a  result 
that  has  since  been  expanded  into  a  generalized  Eliassen- 
Palm  (E-P)  theorem  for  atmospheric  waves  of  various  types 
(e.g.,  Andrews  et  al.,  1987).  In  general  an  E-P-flux  vector 
governs  wave  activity,  with  pu'w'  being  associated  with  its 
vertical  component,  the  most  relevant  component  for  gravity 
waves.  Drag  on  the  mean  flow  caused  by  wave  dissipation  is 
referred  to  as  GWD  and  is  given  by  the  divergence  of  the 
wave’s  E-P  flux.  For  irrotational  gravity  waves,  this  is  given 


by  —  —  (pt t'w'\  which  is  zero  prior  to  saturation  since  p  u'w ' 
p  oz  '  ’ 


is  constant.  Parametrizations  compute  GWD  by  evaluating 
the  vertical  gradients  of  momentum  flux  that  are  produced  by 
saturated  waves:  “saturation”  is  discussed  further  in  the  fol¬ 
lowing  sub-subsection. 

Thus  gravity  wave  momentum  flux  must  be  specified  accu¬ 
rately  at  its  source.  For  mountain  waves,  this  is  achieved  (as 
discussed  in  Section  a  of  this  Appendix)  by  calculating  the 
surface  pressure  drag  across  the  subgrid-scale  mountains  and 
apportioning  it  to  a  total  vertical  flux  of  horizontal  momen¬ 
tum  of  radiated  waves.  The  momentum  flux  generally  takes 
the  form  kPqIVqI  i/0l/z'2,  where  the  subscript  “0”  denotes  the 
surface  or  near-surface  level,  K  is  a  dimensionless  coefficient 
that  varies  with  definitions,  topographic  shape  and  so  on  (e.g., 
Baines,  1995)  and  h'  is  an  elevation  proportional  to  the  stan¬ 
dard  deviation  of  subgrid-scale  orographic  height.  A  variety 
of  similar-looking  expressions,  which  vary  in  details  (e.g.,  the 
treatment  of  low-level  blocking  and  wave -breaking)  have 
been  derived  in  parametrizations  for  various  large-scale  mod¬ 
els:  Kim  and  Arakawa  (1995)  show  a  table  summarizing  most 
of  the  earlier  schemes. 


2  SATURATION 

While  various  atmospheric  processes  can  damp  waves, 
parametrizations  focus  mostly  on  dissipation  due  to  wave- 
induced  instabilities  that  develop  when  wave  amplitudes 
become  sufficiently  large.  This  is  a  complex  non-linear 
process  that  is  not  fully  understood  (see,  the  review  by  Fritts 
and  Werne,  2000),  but  is  characterized  macroscopically  by 
limiting  “saturation”  amplitudes  for  the  waves,  beyond  which 
instabilities  develop  and  gravity-wave  dissipation  and  drag 
occur  (Fritts,  1984,  1989;  Fritts  and  Alexander,  2003). 

GWD  parametrizations  utilize  these  ideas:  many  schemes 
classify  waves  into  a  “saturated”  regime,  where  dissipation 
processes  limit  the  wave  amplitudes  and  produce  drag,  and  an 
“unsaturated”  regime  within  which  waves  do  not  dissipate. 
The  wide  variety  of  parametrizations  that  have  emerged  also 
reflect  many  of  the  current  uncertainties  in  gravity  wave  sat¬ 


uration  theories,  which  continue  to  be  debated.  In  single-wave 
parametrizations,  saturation  is  usually  associated  with  “wave 
breaking”  (e.g.,  Lindzen,  1981;  Palmer  et  al.,  1986; 
McFarlane,  1987),  a  process  somewhat  analogous  to  surface 
ocean  waves  breaking  on  a  shore.  Rather  than  breaking  down 
totally,  however,  it  is  assumed  that  waves  “saturate”  by  dissi¬ 
pating  only  enough  wave  energy  into  turbulent  layers  to 
return  the  wave  to  marginal  stability  (Fritts,  1984,  1989; 
Dunkerton,  1989;  Palmer  et  al.,  1986;  McFarlane,  1987). 

For  spectral  parametrizations,  saturation  is  also  a  key  con¬ 
cept.  However,  its  connotation  varies  due  to  a  wide  variety  of 
proposed  physical  mechanisms  underlying  dissipation  of  and 
interactions  among  a  spectrum  of  waves  (see  Section  4c). 
These  differences  lead  to  quite  different  algorithms  for 
parametrizing  GWD.  For  example,  in  some  schemes  the  onset 
of  wave-induced  instabilities  marks  total  dissipation  and 
removal  of  certain  waves  (e.g.,  Hines,  1997a;  Alexander  and 
Dunkerton,  1999).  In  others,  waves  dissipate  their  amplitudes 
much  more  gradually  with  altitude  (Weinstock,  1990; 
Medvedev  and  Klaassen,  2000;  Warner  and  McIntyre,  2001). 

3  RELATIONSHIP  TO  INTRINSIC  PHASE  SPEED 
Quantitatively,  the  wave’s  intrinsic  phase  speed  c  (given  in 
2D  by  c  =  c  -  U(z)',  see  Section  c  of  this  Appendix)  is  an 
important  quantity  in  specifying  saturation  amplitudes  and 
GWD  in  parametrizations.  “Convective”  saturation  sets  in 
when  waves  attain  an  amplitude  that  causes  potential  temper¬ 
ature  surfaces  to  “overturn”:  d/dz(8  +  0')  <  0,  where  0  is  the 
mean  potential  temperature  and  0'  is  the  gravity-wave  poten¬ 
tial  temperature  perturbation.  For  a  2D  hydrostatic  irrotation¬ 
al  gravity  wave,  this  occurs  when  the  peak  horizontal  velocity 
amplitude  of  the  wave,  upEAK,  exceeds  Icl,  and  so  the  satura¬ 
tion  amplitude  equals  the  intrinsic  phase  speed:  (a'PEAK)SAT  = 
Icl  (Fritts,  1984).  In  more  complex  saturation  schemes,  wave 
saturation  amplitudes  are  often  directly  proportional  to  Icl  also 
(e.g.,  Medvedev  and  Klaassen,  2000).  Thus,  saturation 
becomes  more  likely  in  regions  where  Icl  is  small,  since  waves 
saturate  at  smaller  wave  amplitudes  and  so  become  unstable 
more  easily.  The  strength  of  the  ensuing  GWD  is  also  sensi¬ 
tively  dependent  on  c:  for  example,  in  some  schemes  the 
GWD  is  proportional  to  Icl3  (Lindzen,  1981;  Medvedev  and 
Klaassen,  1995;  McLandress,  1998;  Alexander  and 
Dunkerton,  1999).  For  single-wave  saturation,  this  arises 
because  (u PEAK) sat  =  ^  implies  (w'peAk^sAt=  ±  ^  ^ 2 ^  fr°m 
basic  polarization  relations  for  gravity  waves  (see,  Eq.  (36)  of 
Fritts  and  Alexander,  2003),  and  thus  (u'w')SAT  °e  lc3l 
(Lindzen,  1981). 

The  sign  of  c  is  also  important.  In  essence,  wave  saturation 
tends  to  “drag”  the  wind  towards  the  ground-based  phase 
speed  c  of  the  dissipating  wave  (Lindzen,  1981),  so  that  ifc  < 
0  [c  <  U(z)]  the  GWD  is  negative  (westward)  whereas  if  c  > 
0  [c  >  U(z)]  the  GWD  is  directed  eastward.  This  explains  the 
origin  of  the  term  “gravity  wave  drag”:  the  first  GWD  param¬ 
etrizations  were  for  stationary  orographic  waves  (c  =  0),  dis¬ 
sipation  of  which  always  drags  the  flow  towards  this  zero 
phase  speed  and  thus  always  decelerates  winds. 
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Throughout  this  paper  we  use  the  familiar  term  “gravity 
wave  drag”,  with  the  understanding  that  this  “drag”  has  a 
broader  definition  that  allows  both  acceleration  and  decelera¬ 
tion  (including  reversal)  of  atmospheric  winds  due  to  dissipa¬ 
tion  of  both  stationary  and  non- stationary  gravity  waves. 

4  RELATIONSHIP  TO  CRITICAL  LEVELS 

Since  Icl  — >  0  on  approaching  a  critical  level,  the  previous 
arguments  show  that  saturation  becomes  essentially 
inevitable  near  2D  critical  levels,  since  (u  peak)  sat  ^  ^  ~ >  0 
as  well.  This  2D  saturation  condition  ( u'peak)sat=  ^  can 
re-expressed  as  c  -  [f/(z)  +  u'(x,z,t)]  =0,  where  u'(x,z,t)  is  the 
horizontal  velocity  oscillation  of  the  gravity  wave.  By  analo¬ 
gy  (from  Section  c  of  this  Appendix)  with  the  regular  critical- 
level  criterion  c  -  U(z)  =  0,  we  see  that  convective  saturation 
can  be  viewed  as  the  onset  of  a  “wave-induced  critical  level” 
caused  by  the  large  velocity  amplitudes  u'(x,z,t)  of  the  wave 
itself.  This  concept  is  revisited  in  the  parametrization  of  low- 
level  orographic  GWD  (Section  3a).  Since  c  =  0  for  oro¬ 
graphic  waves,  the  condition  is  that  the  total  flow 
velocity  U(z)  +  u{x,z,t)  vanishes.  Because  of  this,  overturn¬ 
ing  regions  for  low-level  mountain  waves  in  downstream 
regions  (see  Fig.  4)  are  often  associated  with  so-called  “stag¬ 
nation  points”  in  low-level  orographic  flow  in  upstream 
regions  (Smith,  2001). 

For  multi-wave  fields,  this  condition  can  be  generalized 
still  further  to  c  -  [f/(z)  +  U'(x,z,t )  +  u{x,z,t)\  =  0  where 
U'(x,z,t)  denotes  the  horizontal  velocity  oscillations  of  all  the 
other  gravity  waves  in  the  spectrum.  This  formulation  is  key 
to  the  spectral  GWD  parametrization  of  Hines  (1997a, 
1997b),  in  which  gravity  waves  are  saturated  and  eventually 
removed  (“obliterated”)  by  the  “Doppler-spreading”  effects 
of  U'(x,z,t)  fields  due  to  these  other  waves  in  the  spectrum 
acting  on  wave  intrinsic  frequencies  and  vertical  wavenum¬ 
bers  (see  Section  4). 

Appendix  B:  Abbreviations  and  symbols 


Table  1 .  Acronyms/abbreviations 

Abbreviation 

Term 

2D 

two  dimensions/two-dimensional 

3D 

three  dimensions/three-dimensional 

CCM 

Community  Climate  Model 

E-P 

Eliassen-Palm 

GCI 

global  cloud  imagery 

GCM 

general  circulation  model 

GFDL 

Geophysical  Fluid  Dynamics  Laboratory 

GWD 

gravity  wave  drag 

MF  radar 

medium  frequency  radar 

NWP 

numerical  weather  prediction 

QBO 

quasi-biennial  oscillation 

SAO 

semiannual  oscillation 

Table  2.  Mathematical  symbols 


Symbol 

Meaning 

c 

ground-based  zonal/horizontal  phase  speed 

c 

intrinsic  zonal/horizontal  phase  speed 

f 

Coriolis  parameter 

Frd 

(inverse)  Froude  number 

h 

subgrid- scale  orographic  height 

hd 

vertical  displacement  height 

K 

tuning  coefficient 

k 

zonal/horizontal  wavenumber 

e 

Scorer  parameter 

K 

zonal/horizontal  wavelength 

K 

vertical  wavelength 

m 

vertical  wavenumber 

N 

Brunt- Vaisala  frequency 

n 

total  spherical  harmonical  wavenumber 

P 

atmospheric  pressure 

pju 

zonal/horizontal  wavenumber  spectrum  of  quantity  X 

p 

atmospheric  density 

Ri 

Richardson  number 

Rim 

minimum  (wave-modified)  Richardson  number 

X 

vertical  flux  of  zonal/horizontal  momentum 

t 

time 

tpROP 

gravity  wave  propagation  time 

At 

model  time  step 

U 

zonal/horizontal  wind  magnitude 

u 

background  zonal/horizontal  wind  profile 

u' 

zonal/horizontal  velocity  perturbation  of  a  field  of  waves 

u 

zonal/horizontal  velocity  perturbation  of  a  wave 

11  PEAK 

peak  zonal/horizontal  velocity  amplitude  of  a  wave 

1“  PEAK>SAT 

saturated  peak  zonal/horizontal  velocity  amplitude  of  a  wave 

0) 

intrinsic  wave  frequency 

of 

vertical  pressure  velocity  perturbation 

W 

vertical  velocity  perturbation 

w 

Fourier-transformed  vertical  velocity 

X 

zonal/horizontal  displacement 

Ax 

zonal/horizontal  model  grid  size 

z 

altitude 

zc 

critical  level 

A  Z 

vertical  model  grid  size 
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